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ABSTRACT 

The process of tidal dissipation inside Jupiter is not yet understood. Its tidal quality factor (Q) is 
inferred to lie between 10 5 and 10 6 . Having studied the structure and properties of inertial-modes in 
a neutrally buoyant, core-less, uniformly rotating sphere (Wu 2004), we examine here their effects on 
tidal dissipation. The rate of dissipation caused by resonantly excited inertial-modes depends on the 
following three parameters: how well they are coupled to the tidal potential, how strongly they are 
dissipated (by the turbulent viscosity), and how densely distributed they are in frequency. We find 
that as a function of tidal frequency, the Q value exhibits large fluctuations, with its maximum value 
set by the group of inertial-modes that satisfy Su> ~ 7, where Su> is the group's typical offset from an 
exact resonance, and 7 their turbulent damping rates. These are intermediate order inertial-modes 
with wave-number A ~ 60 and they are excited to a small surface displacement amplitude of order 
10 3 cm. The Q value drops much below the maximum value whenever a lower order mode happens 
to be in resonance. In our model, inertial-modes shed their tidally acquired energy very close to the 
surface within a narrow latitudinal zone (the 'singularity belt'), and the tidal luminosity escapes freely 
out of the planet. 

Strength of coupling between the tidal potential and inertial-modes is sensitive to the presence of 
density discontinuities inside Jupiter. In the case of a discreet density jump, as may be caused by the 
transition between metallic and molecular hydrogen, we find a time-averaged Q ~ 10 7 , with a small 
but non-negligible chance (~ 10%) that the current Q value falls within the empirically determined 
range. Whereas when such a jump does not exist, Q ~ 10 9 . Even though it remains unclear whether 
tidal dissipation due to resonant inertial-modes is the correct answer to the problem, it is impressive 
that our simple treatment here already leads to three to five orders of magnitude stronger damping 
than that from the equilibrium tide. 

Moreover, our conclusions are not affected by the presence of a small solid core, a different pre- 
scription for the turbulent viscosity, or nonlinear mode coupling, but they depend critically on the 
static stability in the upper atmosphere of Jupiter. This is currently uncertain. Lastly, we compare 
our results with those from a competing work by Ogilvie & Lin (2004) and discuss the prospect of 
extending this theory to exo-jupiters, which appear to possess Q values similar to that of Jupiter. 
Subject headings: hydrodynamics — waves — planets and satellites: individual (Jupiter) — stars: 
oscillations — stars: rotation — turbulence 



1. INTRODUCTION 
1.1. the Puzzle 

We tackle the classical problem of tidal dissipation in 
Jupiter. In the following, we briefly review the problem, 
both for Jupiter and for close-in extra-solar planets. For 
a contemporary and expansive overview of this issue, in- 
cluding a detailed discussion of previous work, we refer 
the readers to Ogilvie & Lin (2004, hereafter OL). 

As Jupiter spins faster than the orbital motion of its 
nearest satellite (Io), Io raises a time-dependent tide on 
Jupiter. The dissipation of this tide in Jupiter transfers 
its angular momentum to Io and spins down Jupiter. We 
adopt the convention of quantifying the inefficiency of the 
dissipation 1 by a dimensionless quality factor Q, which 
is the ratio between the energy in the (equilibrium) tide 

1 This assumes Q is independent of the orbital phase. Hut (1981) 
and others have adopted instead a constant lag time t = e/(Q — 
/), where Q and / are the rotational and instantaneous orbital 
angular velocity, respectively, and / is the free anomaly. These 
two approaches are comparable if Q is frequency independent and 
if the eccentricity is not too large. 
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where e is called the lag angle. It corresponds to the 
angle between the directions of Io and the tidal bulge 
when we are concerned with the equilibrium tide. The 
rate of tidal synchronization scales inversely linearly with 
Q (Murray & Dermott 1999), 
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where Q is Jupiter's spin frequency, u)\ Q — (GMj/a 3 ) 1 / 2 
is Io's orbital frequency, a its orbital separation, Mj, Rj, 
fta, a are Jupiter's mass, radius, tidal love number (&2 ~ 
0.38) and moment of inertia constant (a = I/MjRj ~ 
0.25), respectively. M\ Q is Io's mass. 

Based on the current resonant configuration of the 
Galilean satellites, Jupiter's Q value has been estimated 
to be 10 5 < Q < 2 x 10 6 with the actual value likely 
closer to the lower limit (Goldreich & Soter 1966; Peale 
& Greenberg 1980). The interior of Jupiter is comprised 
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of (at most) a small heavy-element core, a metallic hy- 
drogen region and a molecular hydrogen envelope (see, 
e.g. Guillot et al. 2004), with convection being the dom- 
inant heat transfer mechanism outside the core. The 
most reliable theoretical estimate for the Q value - based 
on turbulent viscosity acting on the equilibrium tide - 
puts Q <~ 10 13 (Goldreich & Nicholson 1977), well above 
the inferred value. The physical origin for this low Q 
value (and thus higher than expected dissipation) has re- 
mained elusive for a few decades, with suggestions rang- 
ing from a substantial inner core (Dermott 1979), to he- 
lium hysteresis around the depth of hydrogen metallic 
phase transition (Stevenson 1983), to a postulated strati- 
fication in the interior that harbors rotationally-modificd 
gravity- modes (Ioannou & Lindzen 1993). Each proposal 
promises interesting implication for the physics of dense 
matter or for the structure of Jupiter. Where does the 
truth lie? Intriguingly, Saturn's inferred Q value is sim- 
ilar to that of Jupiter (Goldreich & Soter 1966). 

The discovery of close-in extra-solar jupiters has reju- 
venated our interest in this problem and provided new 
insights. While the majority of exo-planets are in ec- 
centric orbits around their host stars, the closest-in ones 
have low or nearly zero eccentricities. This results from 
the dissipation of stellar tide inside the planets which 
converts orbital energy into heat without removing or- 
bital angular momentum. Orbital circularization due to 
tidal dissipation inside the planet proceeds at a rate (Hut 
1981): 
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where k 2 , Q and Rj are the planet's tidal love number, 
tidal quality factor and radius, respectively. It orbits its 
host star (mass M*) with a semi- major axis a and an 
orbital frequency u. Fig. 1 in Wu (2003) shows that 
the observed upper envelope of planet eccentricity as a 
function of its semi-major axis can be explained by a 
tidal quality factor of Q m 3 x 10 5 if these are gaseous 
planets similar to Jupiter in their ages and sizes. 2 

The close-in exo-planets and Jupiter may well have 
different formation history, leading to different core 
sizes and different interior compositions. They certainly 
evolve in very different thermal environments, resulting 
in diverging thermal structure in their upper atmosphere. 
Nevertheless, they share similar Q factors. This prompts 
us to seek a physical explanation for Q which is based on 
overt similarities between these planets. The first trait 
in common which we believe is important is that their in- 
teriors are fully convective. The second trait is that they 
rotate fast. Jupiter spins roughly four times for every 
Io orbit, while the spin of close-in (o < 0.1AU) planets 
should have long been (pseudo-)synchronized with their 
orbital motion. So in both cases, the (dominant) tidal 
forcing frequencies viewed in the planets' rotating frame 
are below 2fi. 3 Could these two common traits be re- 
sponsible for the tidal Q values? 

2 One exception is the planet HD 806066 whose abnormally high 
eccentricity may be acquired relatively recently (Wu & Murray 
2003). 

3 In this respect, it is interesting to point out that tidally cir- 
cularizing solar-type binaries have convective envelopes and likely 
spin fast. Curiously, they exhibit similar Q values as these giant 
planets Mathieu et al. (2004). 



1.2. The Inertial-Mode Approach 

In a spinning and neutrally buoyant fluid sphere, a new 
branch of eigen-modes arise: the inertial-modes. Their 
motion is restored not by pressure or buoyancy, but by 
Coriolis force. In the rotating frame, these modes have 
frequencies ranging from zero to twice the spin frequency. 
As noted above, the tidal frequencies also fall in this 
range. How does the presence of these modes affect tidal 
dissipation? 

We have previously studied inertial-modes in non- 
uniform density spheres (Wu 2004, hereafter Paper I), 
focusing on properties relevant to tidal dissipation. We 
found that inertial-modes which can couple to the tidal 
potential are much denser in frequency space compared 
to gravity- or pressure-modes, allowing for good reso- 
nance with the tidal forcing. Inertial-modes have unique 
"singularity belts" near the surface where both mode am- 
plitudes and velocity shear are the largest, leading to 
strong turbulent dissipation. Both these facts suggest 
that inertial-modes are good candidates to explain the 
tidal dissipation in planets. In this paper, we explore 
this possibility for Jupiter. 

Because of mathematical difficulties, rotation has been 
largely ignored in tidal theories (for an exception, sec 
Savonije et al. 1995, as well as their subsequent papers). 
However, this can not be justified when rotational fre- 
quency is comparable to or faster than the tidal fre- 
quency. Tidal response of the fluid is strongly influenced 
by rotation. Our results here show that when rotation 
is taken into account, even the most rudimentary treat- 
ment gives orders of magnitude stronger tidal dissipation 
than when it is not. 

In this direction, most noteworthy is a recent indepen- 
dent work by OL, which appeared while we were writing 
up our results. In this paper, OL calculated the effect of 
inertial-modes in planets, based on essentially the same 
physical picture as we consider here. We discuss their 
work in the context of our results. For un-initiated read- 
ers, we recommend their excellent and helpful review for 
issues related to tidal dissipation and to inertial-modes. 

1.3. Organization 

Paper I has laid a foundation by studying properties 
of inertial-modes. In §2, we first summarize results from 
that paper, then proceed to discuss two issues of impor- 
tance, i.e., how strongly an inertial-mode is coupled to 
the tidal potential, and how strongly an inertial-mode 
is damped by turbulent viscosity. Relevant contents of 
a highly technical nature are presented in Appendix D, 
where a simple toy model is constructed to help explain 
the results. In §3, we discuss the effects of inertial- 
mode dynamical tide on the tidal Q factor, using equi- 
librium tide as a comparison to illustrate the advantage 
of inertial-modes. Lastly, we discuss uncertainties in our 
model, and compare our results with previous work (§4). 
We summarize and discuss other possible applications in 
§5- 

2. INERTIAL-MODES - RELEVANT PROPERTIES 

In Paper I, we show that by introducing the ellipsoidal 
coordinates (Bryan 1889), the partial differential equa- 
tion governing inertial-modes can be separated into two 
ordinary differential equations, both when the density is 
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uniform (Bryan 1889), and when the density satisfies a 
power-law (p oc [1 — (r/R) 2 } 13 , where R is the planet ra- 
dius and r the spherical radius). Moreover, for spheres 
of smooth but arbitrary density laws, we find that one 
could obtain sufficiently accurate (albeit approximate) 
eigenfunctions using these coordinates. 

Each inertial-mode in a sphere can then be character- 
ized by three quantum numbers, n\,n 2 and m. Here, n\ 
and n 2 are the number of nodes along the X\ or x 2 el- 
lipsoidal coordinate lines, and m is the usual azimuthal 
number. All perturbations satisfy the form e lm ^ with <j> 
being the azimuthal angle. For a graphical presentation 
of an inertial-mode, see Figs. 4 & 5 in Paper I. We also 
introduce in Paper I the dimensionless mode wavenum- 
ber A <~ 2(ni + n 2 ), which is related to the dimensionless 
mode frequency fi = lu/20, ss sin(ni7r/A), where u) is the 
inertial-mode frequency viewed in the rotating frame, Q 
the spin frequency, and < p, < 1. Under this con- 
vention, m < denotes retrograde modes, while m > 
prograde ones. 

2.1. Goodness of Resonance 

In a non-rotating star, each eigenmode is identified by 
three quantum numbers n, £, m where n is the number of 
nodes in the radial direction, and £, m relate to a single 
spherical harmonic function Pp(6,4>) that describes the 
angular dependence of the mode. In contrast, the an- 
gular dependence of each inertial-mode is composed of a 
series of such spherical harmonic functions. This has the 
consequence that while only the £ = 2, m = — 2 branch of 
non-rotating modes can be driven by a potential forcing 
of the form P 2 ~ 2 (the dominant tidal forcing term), ev- 
ery even-parity inertial-mode can potentially be driven. 
In this sense, the frequency spectrum of inertial-modcs 
is dense, and the probability of finding a good frequency 
match (mode frequency w forcing frequency) is much im- 
proved over the non-rotating case. 

For a given forcing frequency po, how far in frequency 
does the closest inertial-mode lie? We limit ourselves to 
inertial-modes with A ~ 2(n\+n 2 ) < A max - Approximate 
mode frequency by p s=s sin(ni7r/A) ~ niir/X. Modes 
with the same n 2 but different ni are spaced in frequency 
by ~ 7r/A > 7r/A max . Now allow n 2 to vary between 1 and 
~ A max /4, 4 we find that the best frequency off-resonance 
to no is typically 
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For comparison, gravity- or pressure-modes in non- 
rotating bodies can at best have a frequency detuning 
of Suj/ui <~ 1/n with n being the radial order for the 
mode of concern. 

2.2. Overlap with Tidal Potential 

Io orbits Jupiter in the equatorial plane with a fre- 
quency uo\ = 2ir / 1.769 day -1 and at a distance a, while 
Jupiter spins with a frequency = 27r/0.413 day - . 
Viewed in Jupiter's rotating frame, Io rotates retro- 
gradely with frequency u' = £1 — lo\ and exerts a pe- 
riodic tidal forcing on Jupiter. We ignore Io's orbital 
eccentricity (e = 0.004) in this problem. So at a point 



(r, 9, 4>) inside Jupiter, the potential of the tidal pertur- 
bation can be decomposed as 
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The first term is necessary for maintaining the Keplerian 
motion of this point in Jupiter; the second term corre- 
sponds to the potential when Io is smeared into a ring 
along its orbit; the third term is the one of relevance 
here. It describes the periodic forcing by Io in Jupiter's 
rotating frame. Keeping only this term and writing 
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we obtain /Li t ide = <^tidc/2fi = 2u>'/2f2 = 0.766 and m = 
—2. Here, w = rsin6> is the cylindrical radius. 

We investigate here the coupling between inertial- 
modes and the above tidal potential. Assuming the two 
have the same time-dependence, we integrate the forcing 
over the planet and over a period to yield the overlap 
work integral, 

j> dt J d 3 r p^ ■ V<5$ tidc = - J d 3 r 6<5> tidc V ■ (p£) 
= - J d 3 r S^ tidcP ' = - J d 3 r6$ 
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Here, £ and p' are the displacement and Eulerian den- 
sity perturbation from the inertial-mode, while its wave- 
function V is related to p' by p 1 = oj 2 p/c 2 ip (eq. [9] 
in Paper I). The overlap integral represents the energy 
pumped into the mode per period. 

2.2.1. Tidal Overlap for the Equilibrium Tide 

In the limit where the tidal frequency falls well be- 
low the dynamical frequency of the planet, the latter 
reacts almost instantaneously to satisfy hydrostatic equi- 
librium. This tidal response is termed the 'equilibrium 
tide'. An extra response arises when Wti de has a near- 
resonant match with one of the free modes in the planet, 
and this is called the 'dynamical tide'. Physically speak- 
ing, the 'equilibrium tide' is the sum of all the 'dynamical 
tide' response driven at off-resonance. 

Tidal overlap for the equilibrium tide is the largest 
among all tidal response. Disregard any time derivative 
in the fluid equation of motion, take TV 2 = for the 
neutrally buoyant interior, and assume any perturbation 
to be adiabatic, we use equations in §2.1 of Paper I to 
obtain the following instantaneous response, 
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This is the energy stored in the equilibrium tide and it 
appears in equation (1). Taking M\ Q = 8.93 x 10 25 g, a = 
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4.22 x 10 10 cm, and adopting a Jupiter model from Guillot 
et al. (2004), we find Eo ~ 3 x 10 30 erg. In comparison, 
the current potential energy of Io is <~ 3 x 10 38 erg. So, 
over the history of the solar system, Jupiter could have 
pushed Io outward for a negligible 10~ 7 of its current 
orbit if Q <~ 10 13 , namely, only 10 -13 fraction of Eo is 
dissipated per tidal period 

As a side note, the spatial dependence of the tidal 
potential, as well as that of the equilibrium tide, can 
be expressed in the following form which resembles the 
spatial dependence of an inertial-mode in a uniform- 
density sphere: <5$tidc °t ^ °t 

P^ 2 (x 1 )P^ 2 (x 2 ), here 
x\ and xi are the afore-mcntioncd ellipsoidal coordi- 
nates. In comparison, the lowest-order incrtial-mode 
(m = n 2 = 0, also called a R-mode) has a wavefunc- 
tion ijj cx P 3 _ (xi)P$ (X2). 

2.2.2. Tidal Overlap for Inertial- Modes 

Consider first the tidal coupling of a gravity-mode in a 
solar-type star. Firstly, this mode needs to have a spher- 
ical degree I = 2 and an azimuthal number |m| — 2 to be 
compatible with the tidal potential. Its radial eigenfunc- 
tion oscillates quickly in the WKB region and flattens 
out in the upper evanescent region (the convection zone) . 
Overlap with the (smooth) tidal potential therefore is 
largely contributed by the evanescent region, with the 
contribution from different nodal patches in the WKB 
region canceling out each other. 

The situation is different for an incrtial-mode. Firstly, 
every even- parity, |m| = 2 inertial-mode contains a i = 2 
spherical component that can couple to the tidal poten- 
tial. Moreover, the upper evanescent region of an inertial 
mode is comparable in size to any other nodal patch but 
with much lower density. As such it is not particularly 
important for the tidal overlap. The net tidal overlap is 
the small residue after the cancellation between all re- 
gions. This property makes it difficult to reliably calcu- 
late the overlap integral. In fact, obtaining results in this 
section has been the most difficult part of this project. 
Much attention is paid to ensure the accuracy of nu- 
merical integrations, and to analytically understand the 
numerical results. 

We delegate much of the technical discussions to the 
appendixes. In appendix §B, we evaluate tidal overlap for 
inertial-modes in a uniform-density model. In appendix 
§C, we discuss results for models of a single power-law 
index (0). Lastly, in appendix §D, we present results 
for models with more realistic density profiles, including 
ones from Jupiter models. We substantiate our numerical 
results by studying a simple toy-model where analytical 
results are available. Here, we list relevant conclusions. 

We find that the severity of cancellation rises with in- 
creasing mode order. We quantify this severity by the 
following dimcnsionless number, 
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While C n = 1 for the equilibrium tide, C n decreases with 
rising A (or with rising n where n = rii + n 2 ) with a 
slope that depends on the model. In detail, integration 
of the top integral over the spherical angles always leads 
to a cancellation of order n -1 , while integration over the 



radius suffers a cancellation with a magnitude that de- 
pends on factors like the polytropic index of the model, 
or discontinuities in density or density gradient. 

As is shown in Appendix B, in a uniform-density 
sphere, tidal overlap for all modes is zero because the 
material is incompressible (p' = 0). When wc adopt a 
constant pressure, constant density sphere, we find that 
only the two lowest order even-parity modes couple to 
the tide (Papaloizou & Savonije 1997). 

For models satisfying a single power-law density profile 
(p oc [1 — (r/i?) 2 ]' 3 ), C n ~ l/n 2l3+1 for even-parity modes. 
For instance, = 1 and = 1.8 yield C n ~ 1/n 3 and 
C n ~ 1/n 4,6 , respectively. This expression is obtained 
from a simple toy-model and is supported by integration 
of the actual inertial-mode eigenfunctions (Appendix C) . 

Inside Jupiter, gas pressure satisfies the ideal gas law 
above a radius r/R ~ 0.98, while it is dominated by 
that from strongly interacting molecules below this ra- 
dius (discussed in Appendix A.l). The density profile 
can be roughly fitted by two power-laws with varying 
from a value of 1.8 near the surface to ~ 1 deeper down. 
This changing affects the tidal overlap. Let the tran- 
sition occur over a radius Ar. We find C n ~ 1/n 3 for 
n < R/Ar and C n ~ 1/n 4 - 6 for larger n values. These 
are expected since lower order modes mostly sample the 
= 1 region and are evanescent in the = 1.8 envelope, 
while higher order modes experience the 0—1.8 power- 
law. Realistic Jupiter models presented by Guillot et al. 
(2004) yield Ar/R — 0.02, or Ar — 4 local pressure scale 
heights. 

The tidal overlap is also affected by discontinuities in 
density or density gradient. The former may occur if, 
for instance, the metallic hydrogen phase transition is 
of the first-order, while latter occurs if it is of second- 
order. For a density discontinuity with a fractional value 
Ap/p, C n ~ (Ap/p) n~ x /n oc 1/n 2 , while for a den- 
sity gradient discontinuity of Ap' /p', the overlap integral 
C n ~ (Ap'/p') 1/n 3 oc 1/n 3 . 

So in conclusion, the magnitude of the cancellation in 
the overlap integral depends on the density profile, both 
its overall scaling with depth as well as its interior dis- 
continuities and sharp changes. 

In Appendix C, we show that one can obtain C n by 
substituting the actual inertial-mode eigenfunction with 
a fast-oscillating cosine function with the same number 
of nodes (see Fig. D10). It is as if one can almost 
make do without detailed knowledge of the eigenfunc- 
tion. This insensitivity leads us to believe that, although 
we are in many cases using an approximate solution for 
the inertial-mode eigenfunction, our results for the over- 
lap integral is reliable (more discussion in Appendix D). 

Why is it necessary to go through all these detailed 
analysis? In the expression for C n , while the denomi- 
nator is fairly straightforward to obtain through direct 
numerical integration, the severe cancellation suffered by 
the integral in the numerator renders the numerical re- 
sults in many cases untrustworthy. For instance, a 10 -4 
inaccuracy in the Jupiter model presents itself as a small 
(but finite) density jump and affects strongly the value 
of C n at large n. 

2.3. Turbulent Dissipation 

We demonstrated in Paper I that energy of an inertial- 
mode is stored mostly in the form of kinetic energy. An 
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inertial-mode causes little compression. As such, its dis- 
sipation is dominated by shear viscosity. 

The viscous force, F„, appears in the equation of mo- 
tion as 

p'$, + 2pVL x £ = - Vp' + ^ p' - pV<5$ + F v , (11) 

P 

where 



where the tidal potential (5$ 
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and 



F v = V 
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and v is the shear viscosity coefficient and arises from 
turbulent convection. We adopt the following mixing- 
lcngth-formula 



^cv^c 



1 + (wt cv /2tt) s ■ 



(13) 



Here v cv , £ cv and t cv are the characteristic convection 
velocity, scale length and turn-over time. When convec- 
tive turn-over time is long relative to the tidal period 
(lot cv 3> 1), the effective viscosity is reduced and we 
adopt a reduction coefficient s to describe this behav- 
ior. We adopt s = 2 in our main study (see Appendix 
A. 2 for more discussion) and discuss in §4 the effects on 
our results when taking s = 1. We further define the 
depth (R — r) at which lut cv /2tt — 1 to be z C rit- For the 
Jupiter models we adopt (see Appendix A. 2 & Fig. A8), 
z crit « 1CT 2 - S i? and 

v-^z/R)- 1 - for z > z crit , 

~2xl0 10 (z/i?) 1 ^ /3 forz<z crit . (14) 

Here, (3 is taken to be the surface value {j3 — 1.8). The 
deeper region where f3 = 1.0 has too weak a viscosity to 
be of concern. 

We assume here that the viscous forcing is small com- 
pared to the restoring force for inertial-modes so we can 
ignore its effect on the structure of inertial-modes. 5 Vis- 
cosity does, however, dissipate mode energy. The rate of 
dissipation is 
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where we have integrated by part taking the surface den- 
sity to be zero. Viscosity always damps so 7 < 0. In the 
following, we consider only the magnitude of 7, so we 
re-define 7 = |7|. 

2.3.1. Dissipation Rate for Equilibrium Tide 

The equilibrium tide suffers turbulent dissipation as 
the tidal bulge rotates around the planet. We calculate 
its rate of dissipation here. 

First, we obtain the displacement function (£) for the 
equilibrium tide. We ignore the effect of rotation here. 
The motion is barotropic so £ is irrotational, we can write 
£ = V[/ r Y£ !m (#, </>)], where f r is a function of radius 
alone. The equation of mass conservation, combined with 
equation (8), yields the following equation for f r : 



1 d 

~J> ~oZ I r 'P 



dr 



dr 



5 This assumption is equivalent of requiring that the rate of 
turbulent dissipation 7 falls much below mode frequency w, an 
assumption we later confirm. 



b = -^J-^-GM^/a 3 . We solve for f r with the follow- 
ing boundary conditions: near the center, the asymp- 
totic expansion of the above equation yields f r ocr 2 , so 
df r /dr = 2f r /r; at the surface, the Lagrangian pressure 
perturbation is zero so £ r = df r /dr = —br 2 /g, where g is 
the surface gravitational acceleration. 

Over the whole planet, f r rises roughly as r 2 , with a 
surface tidal height £ r ~ 60 meters (and a comparable 
tangential displacement). Using the expression of v (eq. 
[14]), we obtain a damping rate of 7 equ i ~ 4 x 10 _16 s _1 . 

This damping is distributed over the bulk of the planet, 
with roughly equal contribution coming from each decade 
of depth (but little from above z CI i t where viscosity turns 
over). This rate can also be estimated using 7 ~ v/R 2 
with v taken to be 10 4 cm 2 /s, the value for the effec- 
tive viscosity at the mid-point of logarithmic depth (see 
Fig. A8). Lastly, this corresponds to an effective Ek- 
man number (ratio of period to viscous time-scale) of 
Ek ~ 10~ 13 . 



2.3.2. Dissipation Rate for Inertial- Modes 

Numerically, it is straight-forward to obtain the dissi- 
pation rates for inertial-modes. It is sensitive only to the 
density profile at the envelope, and is hardly affected by 
phase transition or other density discontinuities in the 
interior. In this section, we first derive how the rate of 
turbulent dissipation scales with inertial-mode wavenum- 
bcr (A w 2(m +n 2 )), and then present numerical confir- 
mations for these analytical scalings, using a variety of 
power-law models as well as realistic Jupiter models. 

We use the WKB properties of inertial-modes, dis- 
cussed in §3.1 of Paper I. Inertial-modes can propagate 
between the center and an upper turning point, defined 
in the (xi,X2) ellipsoidal coordinates by x\ — p. ~ 1/A 
or p — \x2 1 ~ 1/A or both. The physical depth of this 
turning point depends on latitude. At 9 <~ cos -1 p 
(or x\ <~ I £2 1 ~ p), it is closest to the surface with 
zi - 2R/(1 - p 2 )/X 2 - i?/A 2 (the 'singularity belt'); 
while at other latitudes, the depth is ~ R/\. Within 
the WKB cavity, the amplitude of inertial-modes rises as 
l/y 7 ^. In the x\ and x 2 coordinates, nodes are spaced 
by <~ (1 — p)/n\ and p/n 2 , respectively, and each nodal 
patch (in total n\ n 2 of these) contributes comparable 
amount to the total mode energy. 

Viscosity works on the gradient of the displacement. 
An inertial-mode propagates with a roughly constant 
wavelength in most its WKB cavity, but its wavelength 
shrinks drastically near or inside the singularity belt 
(both x\ — p and p— \x 2 \ < 1/A). This is where we expect 
the largest dissipation to occur. To order of magnitude, 
|V£| <~ V 2 tp <~ A 2 i/> in the WKB cavity, while within 
the singularity belt, |V£| - V 2 ip ~ ip/Sx - \ 4 ip. We 
first consider modes for which z\ > z cr j t , 6 so v cx z^ 1 ^ 13 
(cq. [14]) in the region of interest. The work integral of 
turbulent dissipation can be estimated as, 



/ 



d 3 r pv V£ : V£ 



6 For Jupiter, this roughly translates to A < 50 since z cr i t ~ 
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^ QVk^KAW. (17) 

Obviously, the viscous integral is dominated by the con- 
tribution from the belt where z ~ z\ <~ R/X 2 , and where 
~ cos _1 /i. Meanwhile, the mode-energy integral is 
dominated by the WKB cavity with each nodal patch 
contributing a comparable amount, 

/ d 3 rp£ -£oc [ dxi f (x\ - x\)dx 2 p\ Vip\ 2 



<A 2 (W0* 



ni_n 2 oc (pip) 2 Zi nin 2 . 



(18) 



A more accurate scaling for the energy integral has been 
established in Paper I (§3.2), yielding this integral to be 
oc n 2 - 7 oc A 3,5 . This latter scaling is applicable in the 
range of A that is of interest to us and is fairly indepen- 
dent of the density profile. Returning to equation (15), 
we obtain 

7 ocA 3 - 5 + 2 ' 3 . (19) 

Now we consider higher order modes for which z\ < 
z cr j t . Most of the damping still arises from near z\, where 
v oc z 1- ' 3 / 3 . We repeat the scaling exercise in equation 
(17) and obtain 

7 oc A 1 ' 5 ^/ 3 . (20) 

In Fig. 1, we present the numerically obtained damp- 
ing rates for power-law models with (3 ranging from 1.0 to 
3.0. Some of these models have double power-law density 
profiles but only the envelope (3 value affects the scaling 
for the damping rates. 7 These numerical results confirm 
our above analytical scalings. 

We have also computed damping rates using realistic 
Jupiter models published by Guillot et al. (2004). These 
models are discussed in Appendix A and have /3 = 1.8 in 
the outer envelope. The numerical results are shown in 
Fig. 2. They follow the scalings derived above and can 
be summarized as, 



7 = 6 x 10 



=3 x 10" 



-13 



A 
7.59 



7.1 



7.59 



for A < 50, 



for A > 50, (21) 



where we have scaled A by 7.59, the value of A for a low 
order inertial-mode (n\ = n 2 = 1). Even this low order 
inertial-mode is rather more strongly damped than the 
equilibrium tide. Mode with A <~ 50 have z\ ~ z C rit ~ 
10~ 28 i?. Moreover, damping rates depend only on A but 
not on (ni,n 2 ) values. 

7 For these double power-law models as well as for realistic 
Jupiter models, the inertial-mode eigenfunctions arc obtained as 
described in Paper I. 



3. TIDAL Q FOR JUPITER 

3.1. Q value by Equilibrium Tide 

For the equilibrium tide, equation (1) (Goldreich & 
Soter 1966) yields Q C qui = Lo/-f cqu i, where u is the tidal 
frequency in the rotating frame (u) — 2u)' = 2(0 — wi ) = 
1.5320), and 7 eq ui is the turbulent damping rate for 
the equilibrium tide as calculated in §2.3.1. Substitut- 
ing with the value 7oqu i ~ 4 x 10~ 16 s _1 , we obtain 
Qequi ~ 10 12 , while Goldreich & Nicholson (1977) pre- 
sented an estimate of Qequi ~5x 10 13 . The discrepancy 
is partially due to the fact that they have adopted an 
effective < v >^ 10 3 while our effective < v >^ 10 4 
(§2.3.1) - the actual viscosity is of course uncertain, eas- 
ily by a factor of 10. Moreover, their estimate is more 
order-of-magnitude in nature. In any case, dissipation of 
the equilibrium tide, as has been argued long and hard, 
can not be responsible for the outward migration of Io 
and other satellites. 

3.2. Q value by Inertial- Modes 

How much stronger dissipation can inertial-modes 
bring about? Compared to the equilibrium tide, inertial- 
modes have the advantage that they can be resonantly 
driven by the tidal forcing as they are dense in the fre- 
quency range of interest (§2.1), and they are damped 
much more strongly than the equilibrium tide (§2.3.2). 
The disadvantage, however, lies in the generally weak 
coupling between an inertial-mode and the tidal poten- 
tial. Can the first two advantages overcome the last dis- 
advantage? Here, we combine results from previous sec- 
tions to calculate the tidal Q caused by inertial-modes. 

3.2.1. Q value by Individual Modes 

We start by calculating the amount of tidal energy 
dissipated via one inertial-mode. The following forced- 
damped oscillator equation describes the interaction be- 
tween an inertial eigen-mode and the tidal forcing, 



pi + mi + P^lt, = -pVS^tide exp*' 



(22) 



where £ is the displacement, and the three terms on the 
left-hand-side represent, respectively, the inertia, the vis- 
cous damping, and the restoring force. The free mode 
will have an eigenfrcquency of uio- The right-hand-side 
is the tidal forcing with frequency u> which we take to be 
lo w Adopting the substitution £ = a£, p' = ap' with 
the tilded quantities normalized as u> 2 /2 J d 3 rp£ ■ £ = 1, 
multiply both sides by ^ and integrating over the planet, 
we obtain the amplitude a 



C 

a= 2 



c 



(23) 



where the tidal coupling C = J d 3 r p'5$> t idc, the fre- 
quency detuning 5ui = u) — ujq ~ (lo 2 — cjq)/2u, and the 
angle e = tan~ 1 (7/2<5cj) (we assign 7 > for damping). 
For the equilibrium tide (Slo = lo), this angle represents 
the lag-angle between the tidal bulge and the tide-raising 
body, 2e w 2tane = j cqui /Lo = 1/Q cqui (eq. [1] & §3.1). 

Energy in the inertial-mode is simply a 2 , and the en- 
ergy dissipated via this mode over one period can be 
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found by 

AE = I dt^ = (f> dt I d 3 rRe[£\Re[pV5$ tide e iu ' t } 



= \a\C <j> dtuj sm(bjt + S) cos(ui) 



\a\Cir sin S 



(24) 



2(46uj 2 + 7 2 ) ' 
The tidal Q is related to the above quantity by eq. (1) 
2ttE 4£ 7 /4(5w 2 



Q 



(25) 



AE luC 2 \ T 
where again Eq is the energy in the equilibrium tide, and 
the factor in the parenthesis describes the effect of being 
off-resonance. This expression can also be derived more 
simply taking AE = 2ir/uijE = 2ir/ojja 2 . 

We call a mode "in resonance" with the tide when- 
ever 2\6lu\ < 7. The Q factor associated with a resonant 
mode, Qres, is proportional to the dissipation rate and 
inversely proportional to the normalized tidal coupling, 



Qres — 



C 2 LU 



15 u 



J ^4>vj 2 cos(20) d*r 



(26) 

where equation (9) is used. Notice here that all de- 
pendences on Io's mass and semi-major axis drop out, 
leaving only the dependences on the tidal frequency and 
Jupiter's internal structure. 

How does (2 ros behave for different inertial-modes? 
Based on our previous discussions, we introduce the fol- 
lowing scalings with Ao being the wavenumber of refer- 
ence, 



7 = 7o 



-n 



-1.25ra c 



(27) 



Here, C n is the severity of cancellation in the tidal cou- 
pling, expressed by equation (10). The factor 1.25 in the 
second scaling is needed to accurately relate n = ni + n 2 
to A in the range of interest. The normalized tidal cou- 
pling C can be related to C n as 



c=c, 



1 '6*m4£\T p\d 3 r 



,1/2 



oc C n 



oc 



(28) 



U//#lk(nm 2 )V2 
c c 

- OC — . 

n A 

Here, we have used the information that the envelope 
of tp scales as in the WKB region, and that every 

nodal patch in the WKB region contributes comparable 
amount of kinetic energy to the total budget. Again 
zi stands for the upper turning point at latitude 9 = 

COS -1 \l. 

These scalings combine to yield the following expres- 
sion for Q ros : 

Qrcs = QolO {—) = Q070 ( — ) , (29) 



Ar 



Ar 



where Qo is a constant that depends on Jupiter's internal 
structure. 

We obtain numerical results using two realistic Jupiter 
models published by Guillot et al. (2004): models B and 
D. They are discussed in detail in Appendix A. Of par- 
ticular relevance is that, while hydrogen metallic phase 
transition is treated as a smooth transition in model B 
(interpolated equation of state), model D has a first-order 
phase transition and the associated density jump occur- 
ring around r/Rm 0.8. As a result, C n ~ 1/n 3 (n c = 3) 
in model B, 8 while C n <~ l/2n 2 (n c = 2) in model D. 
These scalings are derived analytically in Appendix D, 
and tested using a toy-model integration. In Fig. 2, we 
further demonstrate that these scalings indeed apply to 
inertial-modes, albeit with quite a bit of fluctuations. 

From equation (21), we obtain n 7 = 7.1 for low-order 
modes (A < 50). So hq is expected to be 16.6 for model B 
and 14.1 for model D. We present numerically calculated 
Qrcs in Fig. 2 and they confirm these scalings. Moreover, 
<5 ros ranges from 10~ 4 for the lowest order inertial-modes 
to 10 10 for model B (and 10 8 for model D) when A - 50. 

3.2.2. Overall Q Value 

If we consider multiple inertial-modes each causing Qi, 
the total effect is 



Q = 



1 



(30) 



So at any given tidal frequency, Q is dominated by the 
mode that contributes the smallest Qi. Which mode is 
this and what is the resulting Q value? We derive ana- 
lytical scalings here to answer these questions. 

At a given forcing frequency, Q values (eq. [25]) for dif- 
ferent modes depend on A non-monotonically. Typically, 
as A increases, Q first decreases and then rises sharply. 
This is because low-order modes typically are driven off- 
resonance (2|cW| > 7) while one can easily find high order 
modes to be in resonance with the tide. For low-order 
modes, as A increases, the chance for a good resonance 
with the tidal frequency improves. This compensates for 
the fact that tidal coupling weakens with A and 



4cW 2 2567T 2 ^ 2 Qo 



res 

T 



7o A 4 ) 



(31) 



Here tiq — 2n 1 — 4 = 2.5n c — 2 — n 7 < 0. For high- 
order modes that satisfy 2|<5oj| < 7, increasingly weaker 
tidal coupling accounts for the fact that Q rises with A 
as Q s=a Q ros oc X nQ . The lowest Q value is to be found 
around modes that satisfy 2\Soj\ = 2\u; — u\ <~ 7. This 
occurs at 

A_ _ /l67rf7\ 1/K+2) 
Ao ~ V 7QA 2 , J 



(32) 



For Jupiter models, this yields A ~ 60 (also see Fig. 
2). These are the modes that are most relevant for tidal 
dissipation. They give rise to a minimum Q value 



Q ~ Qoio 



(33) 



8 Two factors contribute comparably to this scaling: the sharp 
transition of equation of state near r/R ~ 0.98 and the discontin- 
uous density gradient at the phase transition point. 
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This roughly corresponds to Q <~ 10 10 for model B, and 
Q ~ 10 8 for model D (see more detailed calculation be- 
low). 

In the following, wc confirm and refine the above an- 
alytical results by a numerical model. While wc have 
a reasonably good handle on mode damping and tidal 
coupling, we do not have a perfect Jupiter model nor 
exact inertial-mode solutions to produce exact mode fre- 
quencies. So we could not reproduce exactly the tidal 
response of Jupiter as a function of Io's orbital period. 
Fortunately, this problem can be circumvented. In the 
following exercise, we produce an artificial spectrum of 
incrtial-modes, with frequencies that satisfy the WKB 
dispersion relation fj, = sinin^ / A) with A w 2(ni + n<i) 
(Paper I) . To this frequency we add a small random com- 
ponent of order dfi = /u/(10ni), which is of order 1/10 
the frequency spacing between neighboring modes of the 
same value. This random component is to encap- 
sulate our above ignorance but neither its size nor its 
sign qualitatively affect our conclusion. In this treat- 
ment, although we will not be able to obtain the exact 
tidal response of the planet at each forcing frequency, wc 
can get a reasonable statistical impression. In fact, this 
is the only logical approach warranted by our current 
knowledge of the interior of Jupiter. 

For each inertial-mode, we assign a damping as in equa- 
tion (21), and a <5 rcs as in equation (29) (different for 
models B & D), and calculate, as a function of tidal fre- 
quency, the Q value for individual modes as well as the 
overall Q value. The results are presented in Figs. 3 & 
4 for the two models. One observes that the Q value 
fluctuates wildly as a function of the forcing frequency. 
While there is a ceiling to the overall Q value, there may 
be occasions when resonance with very low-order modes 
occurs, leading to deep valleys with Q reaching values 
as small as 10. The ceiling, on the other hand, is deter- 
mined by A ~ 60 modes which are always in resonance 
at any forcing frequency. Due to their high <5 ros values, 
modes of orders higher than these are not important. 

The results should be interpreted statistically. One can 
infer from them two pieces of information about Jupiter's 
Q value. The first is the average Q value across a certain 
frequency range, and the second the probability of Q 
value falling below 10 6 in this frequency range. Here, 
the value Q = 10 6 is taken to be the rough upper limit 
for the empirically inferred Q value. 

The definition for the word 'average' deserves some 
deliberation. We follow Goodman & Oh (1997) and 
Terquem et al. (1998) in adopting the following average, 



Q 



(34) 



This is equivalent to a time-weighted average since the 
time a system spends in a certain state is inversely pro- 
portional to the torque at that state. Over the evolution- 
ary timescale, the system quickly moves through the deep 
valleys (large torque) and lingers around in the large Q 
region. This is also where we most expect to find Jupiter 
today. 

We find that for fj, € [0.7, 0.8], Q w 1.4 x 10° for model 
B and Q s=s 5.8 x 10 7 for model D, roughly consistent 
with our analytical estimates. Recall that <5 eq ui ~ 10 12 . 
Moreover, at any forcing frequency, the probability that 



Jupiter has Q < 10 6 is 3% in model B and ~ 10% in 
model D. 

4. DISCUSSION 

Throughout our calculation, we have assumed that 
Jupiter is uniformly rotating, neutrally buoyant and 
core-less. We have also assumed that its internal con- 
vection provides a turbulent viscosity which is quantified 
by the mixing length theory and which is reduced with 
an index s = 2 when the convection turn-over time is 
long compared to the tidal period. We obtained inertial- 
mode eigenfunctions for realistic Jupiter models using a 
combination of WKB approximation and exact surface 
solution (Paper I). 

In this section, we discuss the validity of our various 
assumptions, factors that might influence our results, as 
well as implications of our results. 

4.1. Tidal Overlap 

Firstly, a precaution about tidal overlap. We find that 
this is the trickiest part of our work because inertial- 
modes propagate essentially over the whole planet, with 
a small evanescent region very close to the surface. Re- 
gions of positive and negative tidal coupling lay side by 
side, leading to strong cancellation and extreme sensitiv- 
ity to numerical accuracy. In fact, for a sphere with a 
density profile that follows a single power-law, the net 
tidal coupling decreases with increasing mode order so 
strongly (Appendix C) that numerical precision is soon 
strained even for fairly low-order modes. Inertial-modcs 
are not important for tidal dissipation in these models. 

In a realistic Jupiter model, the cancellation is less 
extreme due to the following two features: the molec- 
ular to metallic hydrogen transition at r/R <~ 0.8 (ei- 
ther a discreet phase transition or a continuous change) 
and the polytropic index change at r/R <~ 0.98 where 
hydrogen molecules change from ideal gas to strongly 
interacting Coulomb gas (discussed in Appendix A.l). 
These two features act as some sort of 'internal reflec- 
tion' for the inertial-modes - their WKB envelopes inside 
and outside of these features differ. This weakens the 
above-mentioned near-perfect cancellation in the overlap 
contribution from different regions and leads to larger 
tidal coupling. This is confirmed by integration using 
both a toy-model (Appendix D) and actual inertial-mode 
eigenfunctions. In this case, tidal dissipation via inertial- 
modes outweighs that due to the equilibrium tide. 

The inertial-mode eigenfunctions for realistic Jupiter 
models are constructed as follows (see also Paper I). We 
first obtain eigenfunctions for a single power-law model 
with the power-law index (/?) determined by that in the 
outer envelope of the Jupiter model. This can be done 
exactly as long as we ignore the Eulerian density per- 
turbation in the equation of motion. 9 We multiply the 
resulting wave-function by a factor y/ p sur {/ p where p su rf 
is the density for the above single power-law and p is the 
actual density. We showed in Paper I that in the WKB 
region, this construction approximates the actual eigen- 
function to order 0(1/ A 2 ), and it is exact in the surface 
evanescent region. 

9 This term is small and its removal from the equation of motion, 
as we discussed in Paper I, does not preclude tidal forcing between 
the tide and the incrtial-modes. 
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The effects of such a non-exact formulation on mode 
cigenfrcquencies will not significantly alter our results 
and its effects on the damping rates are negligible. But 
does it affect our results on tidal coupling, which, as we 
have shown, depends sensitively even on numerical accu- 
racy? A definitive answer may have to come from high 
resolution numerical calculations. But our toy-model 
gives us some confidence that our approach has captured 
the essence of the problem and that our overlap result is 
qualitatively correct. 

4.2. Turbulent Viscosity 

The next issue concerns the turbulent viscosity. We 
have presented the detailed viscosity profile in Appendix 
A. 2 & Fig. A8.. This is calculated based on the mixing- 
length theory which is order-of- magnitude in nature (also 
see eq. [13]). How much does the Q value change 
when the viscosity is raised (or decreased) by a fac- 
tor of, say, 100? The scaling in equation (33) yields 
Q cx 7o 1 -"Q/("^+ 2 ) j otQoc 7 -°- 8 for model B and cx 7( 7°- 5 
for model D. So even a factor of 100 change in the vis- 
cosity causes a change in the Q value that is comparable 
to our numerical accuracy and is not significant. 

Zahn (1977) has advocated a less drastic reduction 
of the turbulent viscosity when the convective turn-over 
time is much longer than the tidal period: s = 1 in equa- 
tion (13). This produces two differences to our results. 
First, the equilibrium tidal Q is reduced to ~ 10 9 as the 
effective viscosity is increased over the bulk of the planet 
by a factor of ~ 10 3 . Inertial-modes also in general expe- 
rience stronger dissipation, with the change more striking 
for low-order modes. Moreover, modes of lower order can 
now satisfy the resonance condition (2|<W| ~ 7) and they 
are the dominant modes for tidal dissipation. However, 
the enhanced 7 also means every mode now has a larger 
(5 ros , as a result, the overall Q factor by inertial-modes 
is hardly modified from that in the s = 2 case (see Fig. 
5). 

4.3. Density Discontinuities 

As our results in Figs. 3 & 4 show, when there exists 
a discreet density jump inside Jupiter, the overall Q fac- 
tor is ~ 10 7 , or ~ 10 2 times smaller than the case when 
there is no jump, with ~ 10% chance that the current 
Q value falls between 10 5 and 10 6 (the empirically in- 
ferred Q range for Jupiter). This dependence on density 
discontinuity deserves explanation. 

It results from a difference in the overlap integral. In 
the jump case, cancellation in the overlap contribution 
coming from different parts of the planet is less severe 
(C n cx 1/n 2 ), while it is more complete in the no-jump 
case (C n oc 1/n 3 ), as is explained using a toy-model in 
Appendix D. In the no-jump case, the C n cx 1/n 3 scaling 
may arise from two causes: a discontinuity in the den- 
sity gradient due to, for instance, a second-order phase 
transition, and a sharp transition in the power law in- 
dex [3 (equivalently, the polytropic index Ti) when the 
equation of state changes. In Jupiter models, the latter 
occurs at r/R <~ 0.98, spanning a range of Ar/R ~ 0.02, 
or ~ 4 local pressure scale heights (Appendix A.l). The 
overall Q factor is little affected if either transition region 
is shifted upward or downward by a few pressure scale 
heights. However, if the second-order phase transition 



does not exist, and if the polytropic transition occurring 
over a range Ar/R ^> 0.02, we expect C n cx 1/n 4 ' 6 and 
the overall Q factor to be much larger. 

Does Jupiter harbor a density jump? 

One possibility is the so-called metallic hydrogen phase 
transition. Our knowledge of the equation of state for 
hydrogen at Mbar level is currently limited. We do not 
know whether the transition from a molecular fluid to a 
conductive fluid (metallic hydrogen) is a plasma phase 
transition (PPT) with a discreet density jump, or a con- 
tinuous process with only a jump in the density gradient. 
And in the case of PPT, we do not know whether the 
actual Jovian adiabat falls below or above the critical 
temperature for a first-order transition (Stevenson, pri- 
vate communication). Plighted by these uncertainties, 
planet modelers have typically chosen to insert (or not 
to insert) by hand a small density jump at the suspected 
PPT location, and then interpolated between very low 
and very high pressures (where we know the equation 
of state well), under certain assumptions, to obtain the 
pressure-density curves around this point. We build our 
analysis on two examples of such models (model B with 
a smooth transition and model D with a jump). Inter- 
estingly, Guillot et al. (2004) showed that among models 
that match all observational constraints on Jupiter, the 
ones with PPT equation of state have larger core mass 
and lower total mass of heavy elements, while the ones 
with smooth interpolated equation of state tend to the 
opposite. 

Another possibility may follow from helium/hydrogen 
phase separation. Whenever the Jovian adiabat falls be- 
low the critical temperature curve for helium immiscibil- 
ity, helium separates from hydrogen and forms helium- 
rich droplets that fall toward the center (Salpeter 1973). 
Due to its cooler interior, this process has proceeded fur- 
ther in Saturn than in Jupiter. But even in Jupiter there 
may be a density jump, or at worst, a jump in density 
gradient, associated with this effect. 

Close-in hot exo-jupiters presumably have higher over- 
all entropy than Jupiter does, as radiation from their 
surface is effectively sealed off by the stellar insulation. 
Their interior temperature is higher at a given pressure. 
Both PPT and helium rain-out are therefore less likely 
to occur in these bodies. 

In summary, current Jupiter models exhibit features 
that warrant Q ~ fO 9 . It is plausible to find a non- 
negligible density jump in the Jovian interior, due either 
to a first-order PPT or helium/hydrogen separation, in 
which case we obtain Q ~ 10 . This, however, is more 
difficult to justify in hot exo-jupiters, compromising our 
initial goal of searching for a universal mechanism. 

4.4. Presence of a Solid Core 

We have assumed here that convection penetrates into 
the center of Jupiter. But it is possible that Jupiter 
does have a solid core. Dermott (1979) pointed out that 
body tide in the (imperfectly elastic) solid core of Jupiter 
with a core quality factor ~ 30 can account for the ob- 
served tidal dissipation. However, this requires a core 
size which is at the upper-end of current determinations 
(r CO ic/R ~ 0.15) as well as a core quality factor which 
is currently unknown. Moreover, the efficiency of such 
a mechanism depends sensitively on the core size and it 
may be unreasonable to expect that exo-jupiters all have 
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core sizes within a narrow range. So we restrict ourselves 
to consider the effect of a core on the tidal Q factor due 
to inertial-modes. 

Incrtial-modes are excluded from the solid-core. For 
an estimate, we retain the inertial-mode eigenfunctions 
calculated for the core-less case, but suppress from the 
core region contribution to mode energy, mode damp- 
ing, and tidal overlap integral. We find no substantial 
difference between this and the core-less case (one can 
also compare results from model B which is core-less and 
model D which has a 10M© core). Contribution from 
the core region to the overlap integral, for instance, is 
insignificant as the radial integrand drops as oc r 6 (eq. 
[C3]): radial dependence of the tidal potential goes as 
r 2 , and inertial-modes are more anelastic (small p') in 
the high density region. 

A more subtle influence of the core, however, may be 
present. While we have been able to separate spatial 
variables and calculate inertial-mode eigenfunctions in 
the ellipsoidal coordinates for core-less models, the pres- 
ence of a spherical core destroys this convenience. The 
inner boundary conditions can no longer be defined along 
constant ellipsoidal coordinate curves and we have to re- 
turn to the original partial differential equations. This is 
analogous to the situation where the Coriolis force breaks 
the symmetry of a spherical star, with the result that the 
angular dependence of an eigen-mode in a rotating star 
can no longer be described by a single spherical harmonic 
but only by a mixture of them. So it is perceivable that, if 
we adopt core-less inertial-mode eigenfunctions as a com- 
plete basis, inertial-mode eigenfunction in the presence 
of a spherical core may be a mixture of these functions. 
This gives us a hint on how to proceed when there is a 
core. It is possible to obtain the mixing ratio and use 
these to calculate new damping rates, mode energy and 
tidal coupling. We conjecture that the mixture becomes 
purer (more dominated by one component) as the core 
size approaches zero. In particular, we expect the mixing 
not to be important when the core size is much smaller 
than a wavelength of the inertial-mode (r core /R <C 1/A). 
We plan to extend our calculation to the solid core case 
in the future. 

The above conjecture seems to be supported by numer- 
ical calculations by Ogilvie (2005). He recovers low-order 
inertial-modes when he decreases the core size. When 
the core size is significant, however, OL's study discov- 
ered something else. Instead of global inertial-modes, 
they found that fluid response to the tidal forcing is con- 
centrated into characteristic rays which become singu- 
larly narrow as viscosity goes to zero. This appears a 
rather different picture from ours and the physical origin 
of these singular rays deserves understanding. 

4.5. Radiative Atmosphere 

We have also assumed that the convection zone ex- 
tends all the way to zero density. This may be unre- 
alistic for Jupiter, and worse still for exo-jupiters. In 
the Jupiter models we adopted, convection gives way to 
radiation just above the photosphere (p = lbar). The re- 
ality is more complicated (also see discussions in Paper 
I). Temperature in the Jovian atmosphere is such that as 
a fluid parcel travels upward, its water content condenses 
and releases latent heat. The resulting adiabatic gradient 
(the 'wet adiabat') depends on the water content and is 



shallower than the one that does not include water con- 
densation (the 'dry adiabat'). So for a given temperature 
profile, a particularly dry parcel can be convectively sta- 
ble. This is consistent with the Galileo probe data which 
indicates stable stratification down to <~ 20bar after en- 
tering a dry spot on Jupiter (Allison & Atkinson 2001). 
Available Jupiter models are at best 1-D representation 
of the 3-D structure, and our results depend critically on 
the temperature structure and turbulent viscosity in the 
upper atmosphere of Jupiter. 

What is the effect of a thin radiative atmosphere on 
inertial-modes? Inertial-modes may not be perfectly re- 
flected near the surface and some of its wave-flux can 
be smuggled out of the convective region in the form 
of gravity-waves. The radiative zone has a peak Brunt- 
Valsala buoyancy frequency 
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which is much higher than the inertial-mode frequencies 
we are interested in (u) <~ 3.5 x 10~ 4 s _1 ). So the relevant 
gravity-wave is high in radial order and is strongly mod- 
ified by rotation, satisfying N ^> oj <~ f2. Such waves can 
be calculated (semi)-analytically under the 'traditional 
approximation' and are called the 'Hough modes'. The 
smuggled wave-flux is subsequently lost in the higher at- 
mosphere where the gravity-wave breaks. This brings 
about enhanced damping to the inertial-mode. Recall 
that the overall Q factor scales roughly as inverse square 
root of the damping rate. So unless the resultant damp- 
ing rate is orders of magnitude above the rate of turbu- 
lent damping, the overall Q factor is little affected. 

There are other ways in which a radiative envelope may 
affect incrtial-modes. The upper-turning point (z/R ~ 
1/A 2 when 9 ~ cos -1 \i and z/R ~ 1/A otherwise) of 
a sufficiently high order inertial-mode may fall near or 
above the convective-radiative interface. When this oc- 
curs, the structure of the inertial-mode is significantly 
modified. The radiative region imposes a different sur- 
face boundary condition on the inertial-mode than the 
one we assume here (vanishing Lagrangian pressure per- 
turbation). This different boundary condition, as is illus- 
trated by the toy model in Appendix D, may give rise to 
much different (likely larger) tidal overlap and therefore 
a different Q (likely smaller) factor (see also §4.8). 

Extra-solar hot jupiters are strongly irradiated by their 
host stars. Their atmosphere is more isothermal lead- 
ing to a substantially thicker radiative envelope (down 
to <~ 30 km below photosphere) than that in Jupiter. 
This envelope may sustain rotationally-modified gravity- 
waves ('Hough Modes') which may be resonantly (if these 
waves are trapped) excited by the tidal potential. It is 
possible that this explains the tidal dissipation in these 
hot jupiters (Lubow et al. 1997). However, inertial- 
modes should still exist and will couple to the tidal po- 
tential even in these planets. The fact that the Q-values 
appear to be similar between the exo-jupiters and our 
Jupiter leads us to suspect that inertial-modes will re- 
main relevant. It is foreseeable, for instance, that these 
planets harbor a new branch of global modes which are 
inertial-mode like in the interior and gravity-mode like 
in the exterior. 

4.6. Where is the tidal energy dissipated? 
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In our picture of resonant inertial-mode tide, most of 
the tidal dissipation occurs very near the surface, where 
both the kinematic viscosity and the velocity shear are 
the largest. In a realistic Jupiter model, the effective tur- 
bulent viscosity peaks at a depth of ~ 60 km (-z cr it, Fig. 
A8 in Appendix A. 2), and decays sharply inward. Mean- 
while, the displacement caused by inertial-modes rises 
outward toward the outer turning point. And the ve- 
locity shear reaches its maximum inside the 'singularity 
belt' (Paper I), which is found to be around 9 w cos -1 fx, 
with an angular extent ~ R/X and a depth <~ R/X 2 . 
For inertial-modes most relevant for tidal dissipation 
(A <~ 60), this depth roughly coincides with the location 
of maximum viscosity. We have confirmed numerically 
that most of the dissipation indeed occur in this shallow 
belt. 

The tidal luminosity in Jupiter is ~ 7 x 
10 20 (10 6 /Q)erg/s. What is the effect of deposit- 
ing this much energy in a shallow layer? We compare 
this against intrinsic Jovian flux of F ~ 5000 erg/ cm 2 / s. 
The total intrinsic luminosity passing through the belt 
is ~ 2irR 2 /XF - 3 x 10 22 crg/s. This is larger than (or 
at worst comparable to) the tidal luminosity. Another 
way of phrasing this is to say that the local thermal 
timescale is shorter than (or at worst comparable to) 
the ratio between local thermal energy and the tidal 
flux. So the belt is expected to be able to get rid of the 
tidal energy without suffering significant modification to 
its structure. 

Angular momentum is also deposited locally. We as- 
sume here that the convection zone is able to diffuse 
the excess angular momentum almost instantaneously 
toward the rest of the planet. However, if convective 
transport is highly anisotropic and prohibits diffusion, 
it is possible that this (negative) angular momentum is 
shored up near the surface and contributes to surface 
meteorology of Jupiter. 

The transiting planet HD209458b is observed to have 
a radius of <~ 1.3i?j (Brown et al. 2001). Its proximity 
to its host star and its currently near-circular orbit raise 
the possibility that its over-size is a result of (past or 
current) tidal dissipation (Gu et al. 2003). However, if 
our theory applies also to these hot jupiters, we would 
expect that the tidal heat is deposited so close to the 
planet surface that it can not be responsible for inflating 
the planet. 10 Moreover, given the short local thermal 
timescale, any change to the planet structure should dis- 
appear once tidal dissipation ceases. 

4.7. Tidal Amplitude and Nonlinearity 

If inertial-modes are resonantly excited to large ampli- 
tudes, they can transfer energy to other inertial-modes in 
the planet and be dissipated by nonlinear mode coupling. 
To see whether this is important, we consider the ampli- 
tude of inertial-modes. This is largest near the surface 
around the 'singularity belt'. When an inertial-mode is 
resonantly excited (|<W| < 7), we obtain a horizontal sur- 
face displacement £ h ~ 10 n (A/7.59)~ 7 cm. While this 
implies extreme amplitudes for low-order modes, they 
only come into resonance rarely. For modes of inter- 

1(1 It is difficult to imagine how entropy deposited near the sur- 
face can be advected inward to raise the entropy level of the entire 
planet. 



est (A ~ 60), the typical surface displacement amplitude 
is ~ 10 3 cm, 11 so the dimcnsionless amplitude (£/Rj) 
is 10~ 7 . Can such an amplitude incur strong nonlinear 
damping? 

At such small amplitudes, nonlinear effects can be well 
described by three-mode couplings. The efficiency of 
this process scales with the amplitudes of the modes 
concerned. The most important nonlinear coupling is 
parametric resonance: when the inertial-mode reaches a 
threshold amplitude, pairs of daughter inertial-modes, at 
half the frequency and with m = — 1, can be paramet- 
rically excited and can grow to significant amplitudes. 
Nonlinear mode coupling then drains energy quickly out 
of the original mode. The threshold dimcnsionless ampli- 
tude is (Landau & Lifshitz 1969; Wu & Goldreich 2001) 
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where k is the coupling coefficient between the parent 
and the daughter pair, 72 the damping rate for the 
daughter modes, and Su> the frequency detuning for this 
resonance. Arras et al. (2003) has studied the coupling 
coefficient for inertial-modes in a uniform density sphere 
and found k < n\j the maximum coupling coefficient 
obtains for daughter pairs that are spatially similar and 
maximally overlap. 12 We adopt their result here. We 
further take Sco = and 72 = 71 ~ 10~ 7 to obtain the 
lowest possible threshold amplitude. For inertial-modes 
of interest, £/i?j| pa ra ~ 10~ 3 . So parametric damping of 
the tidally forced inertial-modes is un-important. 

Another three-mode coupling of consequence is be- 
tween the inertial-mode, itself and a mode at twice 
the frequency (up-conversion) . However, in the case 
of Jupiter, twice the tidal frequency falls outside the 
inertial-mode range. 

Unconsidered here is another form of parametric reso- 
nance: simultaneous excitation of two inertial-modes by 
the tidal potential, with frequencies of the two modes 
summing up to the tidal frequency. We find this to be 
also negligible for Jupiter-Io system, but likely important 
for exo-jupiters. 

4.8. Comparison with Ogilvie & Lin (2004) 

The most relevant work to compare our results against 
is that of OL, which is an independent study that ap- 
peared while we were revising our paper. In their work, 
the same physical picture as that discussed here was con- 
sidered, namely, tidal dissipation in a rotating planet. 
They employed a spectral method to solve the 2-D partial 
differential equations which describe fluid motion forced 
by the tidal potential inside a viscous, anelastic, neutrally 
buoyant, polytropic fluid. This procedure directly yields 
the value of the tidal torque on the planet, without the 
need of a normal mode analysis. The numerical approach 
allows them to include the effect of a solid core, as well 
as that of a radiative envelope. Overall, they concluded 

11 In contrast, the displacement amplitude of the equilibrium 
tide is much larger, £^ ~ 60 m. Q equ i is large, however, because 
the equilibrium tide is dissipated very weakly. 

12 In their normalization, the dimcnsionless amplitude is unity 
when mode energy equals the rotational energy of the sphere. This 
is similar to setting the dimensionless amplitude to be the ratio 
between displacement and radius at the surface. 
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that inertial-waves can provide an efficient mechanism 
for tidal dissipation, and that the tidal Q factor is an er- 
ratic function of the forcing frequency. We concur with 
these major conclusions. 

However, many technical differences exist between the 
two works. To better understand both works, it is illu- 
minating to discuss some of these differences here. 

Firstly, as is mentioned in §4.4, while we obtain global 
inertial-modes which have well defined WKB proper- 
ties and discreet frequencies, OL demonstrated that the 
tidally-forced response of a planet is concentrated into 
characteristic rays which are singular lines in the limit of 
zero viscosity. While viscous dissipation in their case oc- 
curs in regions harboring these rays, our inertial-modes 
are predominately dissipated very near the surface (the 
'singularity belt'). Moreover, although both our Q val- 
ues exhibit large fluctuations as a function of tidal fre- 
quency, the origin of the two may be different - in our 
case, a deep valley indicates a good resonance between 
the tide and a low-order inertial-mode, while the situa- 
tion is less clear in their case. All these differences may 
originate from the presence (absence) of a solid core in 
their (our) study. We are currently investigating the un- 
derlying mathematical explanation for these differences. 
Again, it is interesting to note that as the core size ap- 
proaches zero, inertial-modes seem to reappear (Ogilvie, 
2004, private communication). 

Secondly, OL's results are based on a n = 1 poly- 
trope, for which we find that tidal coupling is vanish- 
ingly small (see Appendix C), 13 and that inertial-modes 
are not important for tidal dissipation. It is currently un- 
clear whether this difference arises from the presence of a 
core or from the presence of a radiative envelope in their 
study. Despite a steep suppression of the tidal overlap 
integrand near the center (integrand cx r 6 ), the presence 
of a solid core may affect tidal overlap in a more sub- 
stantial manner by reflecting inertial-waves and chang- 
ing their mode structure (§4.4). Meanwhile, a surface 
boundary condition specified at a finite density (instead 
of at p = 0) may cause extra tidal coupling (§4.5), as is 
shown by the analysis in Appendix C. This issue is more 
relevant for extra-solar hot jupiters which have deeper 
radiative envelopes. 

Thirdly, OL assumed a constant Ekman number 
throughout the entire planet. Since 



this implies a viscosity v = ujR 2 Ek ~ 2x 10 16 Ek that is 
constant throughout the planet. We have argued that the 
effective viscosity value for the equilibrium tide should 
be of order ~ 10 4 cm 2 / s (§2.3.1), or an effective Ek <~ 
10~ 13 . However, such a weak viscosity is much smaller 
than is currently reachable by a numerical method in a 
reasonable amount of time. Instead, OL have opted for 
an alternative treatment in which they steadily decreased 
the Ekman number from Ek = 10~ 4 to 10~ 7 and argued 
(based both on numerical evidence and on an analytical 
toy-model) that the final Q value is independent of the 
Ekman number. This contrasts with our results that 
Q roughly scales as 7 (§4-2), obtained for realistic 

13 Although we only present results for a /3 = 1 power-law model, 
they apply to a n = 1 polytrope as well since the two behave 
similarly near the surface and near the core. 



viscosity profiles, where 70 is the damping rate for a mode 
of wavenumber Xq. 

To make the comparison more appropriate, we adopt 
a constant viscosity inside the planet and find that mode 
damping rates 7 = 5 x 10~ 9 (i;fc/10- 7 ) (A/7.59) 3 s" 1 in 
the Jupiter model D, while individual mode <2 res value 
also scales linearly with the Ekman number (eq. [26]). 
Applying scalings derived in §3.2.1, we find an over- 
all Q ~ 2.3 x lO^-Eft/lO -7 ) -1 . This value is consis- 
tent with that obtained by OL for Ek ~ 10~ 7 . Mean- 
while, the equilibrium tide gives rise to Qoqui ~ 4x 
10 6 (Ek/10~ 7 )~ 1 . So in models of a constant Ekman 
number, inertial-modes contribute comparably to tidal 
dissipation as does the equilibrium tide, but no better. 
These results are presented in Fig. 6. 

5. SUMMARY 

In a series of two papers (Paper I & this), we have ex- 
amined the physical picture of tidal dissipation via reso- 
nant inertial-modes. This applies to a neutrally-buoyant 
rotating object in which the tidal frequency in the rotat- 
ing frame is less than twice the rotation frequency. 

In Paper I, we first demonstrate that under some 
circumstances (power-law density profiles of the form 
p cx (1— r 2 )' 3 ), the partial differential equations governing 
inertial-modes can be separated into two ordinary differ- 
ential equations with semi-analytical eigenfunctions. We 
also show that this method can be extended to apply to 
more general density profiles, with the price that the so- 
lution is exact in the surface region but only approximate 
in the WKB regime. Nevertheless, this approximate so- 
lution allows us to draw many physical conclusions con- 
cerning inertial-modes, including their spatial character- 
istics, their dispersion relation, their interaction with the 
tidal potential and with turbulent convection. This semi- 
analytical technique gives us an edge over current com- 
putational capabilities, though full confirmation of our 
conclusions may require careful and high-resolution nu- 
merical computation. It is clear from our study that 
any numerical approach would need to be able to resolve 
the so-called "singularity belt" near the surface where 
inertial-modes vary sharply, and that numerical results 
need to be taken cautiously when evaluating the tidal 
overlap. 

In this paper, we discuss the role in tidal dissipation 
played by inertial-modes. This depends on the following 
three parameters: how well coupled an inertial-mode is 
to the tidal potential, how strongly dissipated an inertial- 
mode is by turbulent viscosity, and how densely dis- 
tributed in frequency are the inertial-modes. We have 
obtained all three parameters using both toy models and 
realistic Jupiter models. Low-order inertial-modes, if 
in resonance (Su> < 7, where 8ui is the frequency de- 
tuning between the tidal frequency and the mode fre- 
quency, 7 is the mode damping rate), can dissipate tidal 
energy with Q as small as Q ~ 10. However, such 
a resonance is not guaranteed at all tidal frequencies, 
and the system sweeps through a fortuitously good reso- 
nance with speed. Inertial-modes most relevant for tidal 
dissipation are those satisfying 6v <~ 7, where 6v de- 
creases with mode wave-number as 5lo cx A 2 , and 7 rises 
steeply with mode wave-number. These are inertial- 
modes with wave-numbers A ~ 60 (or total number of 
nodes n — n\ + n 2 ~ 30). At any tidal frequency, one 
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can always find resonance with one such mode. They 
provide the continuum to the Q value, whereas previ- 
ously mentioned good resonances appear as dense valleys 
superposed on this continuum (see, e.g., Fig. 4). 

The continuum Q value depends sensitively on the 
presence of density discontinuities inside Jupiter, as the 
latter influences strongly the magnitude of coupling be- 
tween the tidal potential and inertial-modes. Current 
Jupiter models show a sharp change in the adiabatic in- 
dex near the surface (hydrogen ideal-gas to Coulomb gas 
transition), this warrants a Q value of ~ 10 9 . The pres- 
ence of a discontinuity in density gradient due to a phase 
transition (metallic hydrogen phase transition and/or he- 
lium/hydrogen separation) has the same effect. On the 
other hand, if the phase transition is first-order in na- 
ture and incurs a density jump, Q ~ 10 7 . Our results 
arc uncertain up to perhaps, one order of magnitude. 
But it is already clear that inertial-modes cause much 
stronger dissipation than the equilibrium tide, which 
yields Q equ i ~ 10 12 . In the case of Q ~ 10 7 , there is a 
~ 10% chance that the current Q value falls between 10 5 
and 10 6 (the empirically inferred Q range for Jupiter). 

Our model also builds on the assumption that Jupiter 
is neutrally stratified and turbulent all the way up to the 
photosphere, as turbulent dissipation for inertial-modes 
with A ~ 60 are calculated to arise mostly near or below 
the photospheric scale-height. Effects like water conden- 
sation may alter the static stability in Jupiter's atmo- 
sphere, making the atmospheric stratification a function 
of space and time. 

We also restrict ourselves to core-less Jupiter models. 
Our conclusion is little affected when we include an in- 
ner core with a size that is compatible with current con- 
straints. However, this is assuming that global inertial- 
modes still exist in the presence of a solid core. Ogilvic 
(2005) extended the study in Ogilvie & Lin (2004) and 
demonstrated that a new kind of tidal response appears 
when Jupiter has a core: fluid motion is tightly squeezed 
into 'characteristic rays' which becomes singular when 
the viscosity goes to zero. This is a drastically different 
picture than the global eigenmode picture described here 
and may lead to different Q factors. 

We have adopted the Goldreich & Keeley (1977) pre- 
scription (s = 2) to account for the reduction in turbulent 
viscosity when the convective turn-over time is long rela- 
tive to the forcing period. Calculations adopting Zahn's 
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prescription (s = 1) produce no difference in the Q value 
caused by inertial-modes, though we find the equilibrium 
tide is significantly more strongly damped. Concerning 
possible effects of nonlinearity: The surface movement of 
inertial-modes is predominately horizontal. For inertial- 
modes that are most relevant for tidal dissipation, the 
surface displacement amplitude ~ 10 3 cm, or ~ 10~ 7 of 
the radius. We estimate that nonlinear effects are negli- 
gible. 

In our theory, tidal heat is deposited extremely close to 
the planet surface (inside the 'singularity belt') and can 
be lost quickly to the outside. For Jupiter, the tidal lumi- 
nosity in this region is smaller than (or at worst compa- 
rable to) the intrinsic luminosity and so would not much 
alter the structure. However, there remains the intrigu- 
ing possibility that the negative angular momentum de- 
posited to the belt may affect surface meteorology (jet 
streams and anticyclones). Moreover, if this theory also 
applies to hot exo-jupiters, the tidal luminosity is un- 
likely to be responsible for inflating planets and solving 
the size-problem of close-in exo-jupiter HD209458b. 

Although our investigation was stimulated by the fact 
that exo-solar planets exhibit similar Q values as Jupiter 
does, it may be difficult to draw a close analogy between 
Jupiter and hot exo-jupiters: the existence of a first- 
order phase transition is less convincing in the latter due 
to their hotter interiors; the upper atmosphere of these 
planets are strongly irradiated by their host stars and 
are therefore likely to be radiative; they may have rather 
different core sizes depending on their formation history. 
Nevertheless, it is our plan to extend the current study 
to exo-jupiters, as investigations into these bodies may 
ultimately yield clue for the story of Jupiter. It is also 
foreseeable that the theory developed here has implica- 
tions for Saturn, Uranus, solar-type binaries, M-dwarfs 
and brown-dwarfs. 
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Fig. 1. — Numerically computed turbulent damping rates for inertial-modcs in various power-law models, as a function of the mode 
wavenumbcr A. Here, we have taken the reduction number s = 2 and included damping rates only for modes with ni ~ ri2, though other 
modes satisfy the same scalings observed here. Models a & b are single power-law models, while c & d are double power-law models 
the index in the interior and P2 that in the envelope) with the transition of (i occurring around r/R ~ 0.98. All models are normalized 
to have the same central density and their viscosity profiles are described by equations (A2)-(A3). The dotted lines are power-law fits to 
the numerical results with the numerical scalings summarized in the top-left corner. These are well reproduced by our analytically derived 
relation 7 oc 

A 3.5+2/3 2 ( eq [ 19 j) Notice that only the envelope power-law index enters the relation. Damping rates in all models flatten at 
large A, and scale with A roughly as A 3 (eq. [20]). 
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Fig. 2. — Tidal coupling, viscous damping rate and resonant Q reB factor for various inertial-modes calculated using two Jupiter models. 
Model B (open circles) has no heavy metal core and no first-order metallic hydrogen phase transition, while model D (solid triangles) has 
a core as well as a density jump at r/R ~ 0.8 due to the plasma phase transition. The upper-left panel presents the (normalized) coupling 
integral C n (eq. [10]) as a function of inertial-mode node numbers (n = n\ +112). Although the scatter is large, model B results arc best fit 
by C n ~ 1/n 3 (solid line), while model D results follow C n ~ l/2n 2 (dotted line). The lower- left panel shows the energy damping rate as 
a function of mode wavenumber A (A ~ 2n). Results from both models scale as A 71 for low-order modes and as A 3 for high-order modes 
(two solid lines), consistent with analytical expectations (§2.3.2). The dotted line in the same panel is the minimum frequency detuning 
as a function of A (eq. [4]). Q res , the Q value contributed by each mode when it is in resonance with the tide (eq. [26]), is plotted on 
the right-hand panel as a function of A. Again, analytical expectations for models B & D are depicted by the solid and dotted lines, 
respectively. While low-order modes (A < 40) from the two models largely share similar Q IeB values, higher order modes follow more closely 
the analytical scalings. Here, we have included only inertial-modes with fi ~ 0.776 but the results remain similar for other inertial-modes. 
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Fig. 3. — Numerically calculated values of Q for model B. The upper panel shows Q as a function of the tidal frequency in the rotating 
frame. Deep dips occur whenever the tide is in resonance with a low-order incrtial-modc (A <C 60), and the ceiling to the Q value is 
determined by the group of modes with A ~ 60 which satisfy 2\8lu\ ~ 7. The tide is always in resonance with one of these modes at any 
frequency. The lower left panel shows the cumulative probability distribution of the Q value within the frequency range 0.7 < fi < 0.8. At 
a given frequency there is ~ 3 percent chance that we will find Q < 10 6 . The probability for this to occur at a given instant in time is 
smaller. The dashed vertical curve locates the time-weighted average Q value (Q, eq. [34]). We find Q = 1.4 X lCr within this frequency 
range. The lower right panel expands the view of the upper panel over this frequency range. The locations of the fine structure in this plot 
are not to be taken literally as we have adopted an approximate dispersion relation for the inertial-modes. 




Fig. 4. — Same as Fig. 3, but for model D in which the tidal coupling decreases as l/2n 2 as opposed to 1/n 3 . This is related to the 
presence of a first-order phase transition at r/R ~ 0.8. While showing overall similar characteristics as those in Fig. 3, Q has now been 
reduced to ~ 5.8 X 10 7 between fi = 0.7 and 0.8, and at any given tidal frequency, there is a ~ 10% chance that Jupiter exhibits Q < 10 6 . 
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Fig. 5. — Same as Fig. 4, but calculated for model D when the index for viscosity reduction is taken to be s = 1 (Zahn 1977) instead of 
s = 2. Mode damping rates now behave as 7 = 10 _1 ° (A/7.59) 4 ' 5 s _1 . We obtain Qsj2x 10 7 , with ~ 30% chance that Q < 10 6 for the 
current tidal frequency. Moreover, Q e qui ~ 10 9 in this case. 
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Fig. 6. — We repeat our calculation for model D, taking the turbulent viscosity to be a constant throughout the planet with the Ekman 
number Ek = u/ujR? = 10~~ 7 . Mode damping rates scale much less steeply with inertial-modc wave- numbers, 7 = 5 x 10~ 9 (A/7.59) 3 ' s" 1 . 
The resulting Q value from inertial-modcs is plotted against fj, in the upper panel, with Q 2.3 X 10 6 over the range ft G [0.7,0.8]. A 
similar calculation for model B yields Q Ri 3.5 X 10 6 . We find that the Q value is inversely proportional to the Ekman number. The lower 
panel translates the Q result into a quantity used in Fig. A2 of OL (the dimcnsionlcss viscous dissipation rate, oc 1/Q), plotted here as a 
function of 2fi = w/Q. The two overlaid lines with Q = 10 (dotted) and Q = 4x 10 (dot-dashed), respectively represent the empirically 
inferred Q value for Jupiter and the Q value associated with the equilibrium tide in this model. These results resemble those presented in 
Fig. A2 of OL for the same Ekman number. 
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APPENDIX 

RELEVANT PROPERTIES IN A JUPITER MODEL 

Here, we study properties of Jupiter that are relevant for the tidal process. This is based on publicly available 
models of Jupiter presented in Guillot et al. (2004). They are produced with the newest equation of state and 
opacity calculations, including the effect of hydrogen phase transition, and alkali metal opacity. They satisfy gravity 
measurements (esp. J 2 & J4) to much better than a percent and reproduce other global properties of Jupiter (radius, 
surface temperature, intrinsic flux). 

In this study, we focus on two models, B and D, out of the five sample models presented in Guillot ct al. (2004). 
Model B is produced with an interpolated hydrogen equation of state (meaning no first-order metallic hydrogen phase 
transition), and has no heavy metal core. Model D, in contrast, contains a first-order phase transition (PPT equation 
of state) and has a core with mass lOM®. The photosphere for both these models is located at a radius of « 7 x 10 9 cm, 
at a pressure of « 10 6 dyne/ cm 2 , and with a temperature 170 K and a density 1.6 x 10~ 4 g/ cm 3 . 

The interiors of these models are fully convective (outside the core). Due to the high density in Jupiter (mean 
density ~ 1.3 g/ cm 3 ), the convection speed needed to carry the small intrinsic flux (5.4 x 10 3 erg/s/ cm 2 ) is highly 
subsonic, resulting in an almost exactly adiabatic temperature profile (super-adiabatic gradient ~ 10~ 8 or smaller). 
This justifies our assumption of neutrally buoyant fluid when investigating inertial-modes. Only the thin atmosphere 
above the photosphere, with a local pressure scale height ~ 20 km, is radiative. 

Density Profile 

Two features in the density profile of these models deserve attention. 

At radius r/R <~ 0.8, pressure ~ 10 12 dyne/ cm 2 , and density <~ lg/cm 3 , hydrogen undergoes a phase transition. 
Above this layer, hydrogen is mostly neutral and molecular. Below this layer, the mean atomic spacing becomes 
smaller than a Bohr radius and electrons are pressure ionized. The strong Coulomb interaction and electron degeneracy 
resemble those in a metal and the transition is referred to as 'liquid metallic hydrogen' transition (Guillot et al. 2004). 
The nature of this transition is still poorly understood. Model B assumes this transition is of second-order and entails 
a discontinuity only in the gradient of density (of order 50%) , while model D assumes it is a first-order transition with 
a density jump of order 10%. These two different treatments should bracket the actual equation of state of hydrogen. 

Another feature sets in nearer the surface, at radius r/R ~ 0.98, pressure ~ 10 10 dyne/ cm 2 and density 0.1 g/ cm . 
Above this region, the gas can be considered as ideal diatomic gas (H 2 ). As the temperature is below 2000 K, the mean 
degree of freedom for each molecule is 5 (three translational plus two rotational). 14 The specific heat per molecule at 
constant volume and constant pressure are, respectively, Cy = 5/2fce, C p = 7/2fce, yielding Ti = <91nP/<9m/9| s = 
Cp/Cy = 1-4. Below this region, however, Ti rises to <~ 1.8 — 2.2 in the main body of the planet, and approaches 3 in 
the very deep interior (Stevenson 1978, 1982). 

This results in different density profiles above and below this region. Recall our definition of /?: p oc [1 — (r/R) 2 } 13 . 
The Jupiter models show that (3 ~ 1.8 (corresponding to Ti ~ 1.4) above this layer, while /3 <~ 1 (corresponding to 
Li <~ 2) in the interior. We also observe that this transition of /? occurs over a fairly narrow region of radial extent 
Ar ~ 0.02_R, or <~ 4 local pressure scale height. As is discussed in §2.2.2, this transition is of significance to our tidal 
coupling scenario. 

But what is the cause behind the rise of Ti near p ~ 10 10 dyne/ cm 2 ? The ionization fraction of electron is too low 
(~ 10~ 6 ) in this region to make a difference by degeneracy pressure; hydrogen is bound into H 2 and only starts to be 
dissociated near p ~ 10 12 dyne/ cm 2 . The true cause, it turns out, is the non-ideal behavior of molecules, a little-talked 
about effect. At a density of 0.1 g/ cm 3 , the mean molecular spacing is <~ 2 A. While the interaction potential between 
H 2 and H 2 molecules is mildly attractive at spacing > 3 A(the van der Waals force), it rises exponentially inward. 
By the time the spacing decreases to below ~ 2 A, this potential is more positive than ksT and the gas pressure is 
no longer dominated by thermal pressure, but is dominated by the repulsive interaction between molecules. This is 
illustrated in Fig. A7. As density rises, molecules increasingly resemble hard spheres, leading to a steeper dependence 
of pressure on density, or Ti <~ 2 (/3 <~ 1). This non-ideal effect loses out at p <~ 10 12 dyne/ cm 2 above which H 2 
molecules are dissociated and electrons are pressure ionized (the metallic hydrogen phase), s 

Turbulent Viscosity Profile 

Inside Jupiter, molecular viscosity is too weak to cause any discernible dissipation on the inertial-modes. We turn 
to turbulent viscosity. 

The kinematic shear viscosity is estimated from the mixing length theory as (Goldreich & Keeley 1977; Zahn 1977; 
Terquem et al. 1998) 

vt ~ Wcv4v— — j^-tt, (Al) 

1 + (wTcv/27r) s 

where v cv , £ cv and r cv are characteristic convection velocity, scale length and turn-over time. The exponent s describes 
the reduction in efficiency when convection is slow compared to the tidal period (lot cv 1). Its value is still under 

14 This number is smaller near the photosphere when the temperature cools toward the rotational temperature of H2 (85 K). Not all 
rotational levels are populated (Saumon ct al. 1995). 
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Fig. A7. — The effect of H2-H2 interaction on gas pressure. The upper panel shows pressure as a function of mean molecular separation 
(in A), with lower density to the right. The solid line is the actual Jupiter profile in a model from Guillot et al. (2004). The dotted line 
represents the ideal gas contribution (thermal pressure) - it falls short of explaining the total pressure above a pressure p ~ 10 10 dyne/ cm . 
The dashed curve shows the contribution from electron degeneracy which only becomes important for p > 10 12 dyne/ cm 2 . For the in- 
between region, another pressure contribution has to kick in. The lower plot examines what this extra contribution is. Here, gas temperature 
inside Jupiter is plotted as a function of the mean separation (solid curve), while the dashed curve depicts the inter- particle potential in 
unit of Kelvin (the <f>YR potential from Ross et al. 1983). Molecular interaction is repulsive for a separation below ~ 3 Aand the interaction 
energy becomes comparable to the thermal energy at a separation ~ 2 A. This contributes to the gas pressure. As density rises, the 
increasingly repulsive interaction dominates the gas pressure and causes the pressure to rise with density more steeply than that of an ideal 
gas. 
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debate, but simple physical arguments (Goldreich & Nicholson 1977; Goodman & Oh 1997) have suggested that s = 2, 
while Zahn (1977) advocated for a less severe reduction with s = 1. We adopt s = 2 in our main study but discuss the 
scenario when s = 1. Some previous studies have adopted a form without 2ir in the above expression. The viscosity 
is effectively smaller but we will show that this does not affect the final Q-value significantly. 

In mixing length theory, v cv « (flux/p) 1 / 3 , r cv ss £ cv /v cv , and £ cv zj ' (3, where H is the density scale height, 
z is the physical depth (z = R — r), and (3 appears in the density power-law as p = [1 — (r/R) 2 ] 13 cx z 13 for z <C R- Let 
the depth at which u;t cv /2tt m 1 be z cr i t . Above z CT n, vt depends on z weakly, 

vr oc z 1 "' 3 / 3 , (A2) 

while below this layer, the turbulent viscosity is significantly reduced and vt decreases sharply inward as 

v T oc z- 1 -^ (A3) 

when 8 = 2 and 

v T oc z" 2 ^ 3 , (A4) 

when s = 1. These approximate scalings are shown in Fig. A8 for Jupiter models B & D. They compare well with 
numerical results. 

TIDAL OVERLAP IN A CONSTANT DENSITY SPHERE 
In a constant density sphere, m = —2 inertial-modes are expressed in the following form (Paper I) 

p' = ^ = AR 2 ^PI\^) Pt\^\ (Bl) 

where A <C 1 stands for the dimensionless amplitude of i[>, and R is the radius of Jupiter. This density perturbation 
is related to the equilibrium tide value p' equi as 

l^^A(l-u 2 )g 1 (x 1 )g 2 (x 2 ) 1 (B2) 



cqui 



3 < Mio 



where Io's orbital frequency uj\ = (GMj/a 3 ) 1 ^ 2 , Mj and Mi Q are the masses of Jupiter and Io, and the dimensionless 
frequency p = oj/2VL = 0.766. The function gi(xi) = P™(xi)/(\ — xf) (introduced in Paper I). 

Pressure in a constant density (p = po = const), self-gravitating sphere is given by p — po[l — (r/R) 2 ] where 
Pa = 2ir/3 GR 2 pf ) = 3/ (8ir)GM 2 /R 4 with M being the total mass. Since [1 -(r/R) 2 } = (x 2 -p 2 )(p 2 -x 2 )/(l-p 2 )/p 2 , 
and volume elements in Cartesian coordinates and ellipsoidal coordinates are related to each other as dxdydz = 
(x\ — x\)/(l — p 2 )/ pdx\dxid4>, we obtain the following tidal overlap, 

f r /.s 9u 2 R 5 M lo A f f 1 p (1 -x 2 )(l -x\)(x 2 -x\) 2/ x 2 , N , , 

The spatial integration can be symbolically performed by Mathematica (best done after conversion to spherical coor- 
dinates) and it yields ~ 0.4(1 — p 2 )/\ where A 2 = 1(1 + 1) — |m|(|m| + 1). So the overlap is 



/ 



S^meP'd'r = 0.4 - -j—T 1 A ■ ( B4 ) 

However, the constant density case is pathological: the value of Ti formally approaches infinity for incompressible 
fluid. Inertial-modes could not cause any density fluctuation (eq. [Bl]) and the tidal overlap is formally zero. 15 

If we take p = constant over the entire sphere (so Y\ is a finite constant), only two motion have non-zero overlap 
with the tidal potential: the equilibrium tide and the two lowest order even-parity inertial-modes with £ = 4. This fact 
has been pointed out in Papaloizou & Savonije (1997) when they considered the convective core of early-type stars. 

TIDAL OVERLAP IN A SINGLE POWER-LAW MODEL 
Are inertial-modes in power-law models coupled to the tidal potential? 

In paper I, we show that one can obtain exact solutions for inertial-modes when the density profile is a single power- 
law p cx [1 — (r/R 2 )] 13 . This allows us to show that inertial-modes in single power-law models do not couple appreciably 
to the tidal potential, except for the two lowest order even-parity modes (corresponding to the t = 4 modes in the 
constant density case). 16 Moreover, the coupling strength falls off with increasing mode order as a power-law with the 
index related to the polytrope index. 

15 The equilibrium tide, on the other hand, has finite tidal overlap. It is equivalent to an inertial-modc with £ = 2 so its spatial overlap 
diverges near the surface as p approaches 0, counteracting the formally infinite Ti . 

16 If we adopt conventional polytrope models with p oc @ , we can obtain approximate solution for the inertial-modes (Paper I). We 
find that they give essentially the same tidal overlap results as single power-law models of the same f3. 
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Fig. A8. — The effective turbulent viscosity vt is plotted here as fine dots against depth (z) for the Jupiter model, when various s values 
are adopted. The dashed curve (s = 0) is the un-reduced turbulent viscosity (corresponding to s = 0). The reduced viscosity (dotted curves) 
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i R r~j 10' cm at which uit cv /2tt ~ 1. Above 



the viscosity is well described by line 

1-/3 „ -,-2.8 



deviate from this curve below a depth z CT - lt 

A: vt ~ 2 X 10 10 (z/i?) 1_ ' 3/ ' 3 oc z 0A (with /3 = 1.8 in the model). Below this depth, reduction is important and vt ~ 4(2//?)^^^^ cx: z 
for s = 2 (straight line B) and i<t~3x 10 5 (2:/i?) — 2 ' 3 / 3 oc z~ 1,2 for s = 1. Deeper down (z > 10 8 cm), as /3 value is varied from 1.8 to 1 
vt takes on a different scaling with depth. However, this is irrelevant as turbulent dissipation from the deep interior is insignificant. 
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The angular dependence of each even-parity, m — — 2 inertial-mode can be decomposed into 

oo 

^ = Mxi)M*2) = J2 p e~ 2 (0)Ce(r), (CI) 

1=1 

where 

C/(r)= J M^i)M^)Pi 2 (O)sm0d0 (C2) 

is non-zero for £ = 2. Ci{r) is an oscillating function of the radius r. We find numerically that C*2(r) oc r 2 near the 
center, while near the surface C 2 (r) approaches a constant for /3 > 1, and oc [1 — (r/i?)]^ _1 ^^ +1 ^ for < (3 < 1. The 
tidal overlap integral is reduced to the following radial integral, 



f cu 2 p 2 3 /32T 3GM Io 2 f* [^3GM± 2 [ R p 2 . 



(C3) 



where we have introduced the integrand Z^M = r A C r (r)p 2 /T\p. It is also an oscillating function of r with an envelope 
that scales as r 6 near the center, and scales near the surface as [1 — (r/i?)]' 3 ~ 1 for (3 > 1, and as [1 — (r/Z?)]^~ 2// ^ +1 ) 
for < /3 < 1. So this integral diverges near the surface if /? < \/2 — 1. 

We find that the integral decreases with increasing mode order in a power-law fashion with the index depending on 
(3. In the following, we explain the observed fall-off with a simple toy-model. 

We approximate D 2 (r) as a product of a rapidly oscillating function and a slowly varying envelope. A rather accurate 
form turns out to be 

D 2 (r)dr = cos(n9)/(9)d9, (C4) 

where the new variable O = cos -1 r/R, n is an integer and is the number of radial nodes in D2(r). We find n = ri\ + n 2 
for the inertial-modes. The smooth function /(6) has a leading term of (tt/2 — 9) 6 near the center (9 <~ 7r/2) and 
a leading term of 9 2/3_1 near the surface (9 <~ 0). 17 . For the moment we assume 2(3 is an integer, and that terms of 
order <d 2/3 and higher also exist near the surface. 
Integrating-by-part yields 

r R rlt/2 



\ D 2 (r)dr = - / cos(n9)/(9)d9 
Jo Jo 



sin(n9)/(9) 



n 



71,2 cos(n9)/'(9) 



tt/2 



+ 



sin(n9)/"(6) 



tt/2 

+ 





cos(n9)/"'(9) 



tt/2 



1 



So the value of this integral depends only on behavior of the function /(9) at the two boundaries. When n is an 
even integer, only odd-order derivatives enter the above expression and we obtain the following results for the tidal 
integral, 



f 

Jo 



f(20-l+Mod[2P,2]) (q\ 

D2(r)dr~ J n2 , +Modm2] [ ' , if 2/3 < 7 

f(7) (_) 

w- if 2(3 > 7 (C6) 

where /^ 2 ^ _1 -'(0) = d 2/3_1 //ii9 2 ^ _1 |e=o and so on. When 2(3 is odd, the above scaling depends on the fact that near 
the surface, terms scaled as 9 2/3 and higher also exist. If they do not (as in the left panel of Fig. C9), 1/n 8 scaling 
prevails. 

When n is an odd integer, slightly different scalings apply: 

R y(2/3-l+Mod[2/3-l,2]) (n\ 

D2(r)dr* J - n2p+Modm _ h2] [ , if 2/3 < 7 

f(6) (_) 

^- if 2/3 > 7 (C7) 

n' 

We have confirmed these scalings numerically with a range of expressions for /(9). The result only depend on the 
boundary behavior of /(9) as long as it is sufficiently smooth. 18 This explains why models with different polytrope 
representations (p oc [1 — (r/R)] 13 or p oc give rise to essentially the same overlap integrals. Moreover, when 

(3 is a fractional number (other than an integer or a half- integer) , we find numerically that J Q R D 2 (r)dr oc 1/n 2 ^ for 
2/3 < 7, similar to the above expressions. 

17 Here, we focus only on models with /3 > 1 

18 In Appendix D, we discuss what the meaning of 'sufficiently smooth' is. 
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Fig. C9. — Severity of cancellation in the overlap integral as a function of mode nodal number (n where n is even) in three single 
power-law models (solid triangles for (3 = 1.0 solid circles for 1.5 and open circles for 1.8). The left-hand panel is the toy model result 
where we have taken the envelope of the cosine function (eq. [C4]) to be /(©) = / = r p /pdr/dQ where © = cos r/R. This allows 
the toy model tidal integrand to have the correct asymptotic behavior as the realistic tidal integrand both near the center and near the 
surface. The severity of cancellation is measured here by J /cos(ri©)d©/ J f\ cos(n0)|d© and it scales as (solid lines) n~ 2 , nT 8 and n~ 3,6 , 
respectively, for the three models, consistent with results in equation (C6). The right panel is the severity of cancellation C n (eq. [10]) 
calculated for incrtial-modcs in the same three models. Again, the three straight lines are the analytically expected scalings, n , n~° 
and n~ 4 ' 6 , respectively, for the three models. The extra power of n compared to those for the toy model arises from cancellation in the 
angular direction, except for the fi = 1.5 model, which does not fall off as n~ 9 due to the presence of 3 term near the surface. Results in 
the j3 = 1.0 model first deviates from the scaling but returns to it at large n and the fi = 1.5 model falls off more steeply than the fi = 1.8 
model, as is expected. 



Recall that the angular integration to yield C2(V) already involves a cancellation of order 1/n. 19 Moreover, even- 
parity modes implies n — n\ +ri2 to be an even number. So for the following three power-law models, = 1.0, = 1.5 
and — 1.8, we expect that the overall tidal overlap falls off with n as n~ 3 , n~ 5 and n~ 4 - 6 , respectively. These 
analytical expectations are plotted in Fig. C9 along with numerical results. The agreement is reasonable, both when 
integrating using the toy model (/(©) cos(nO)) and when integrating using realistic inertial-mode eigenfunctions. 

In obtaining results like those presented in Fig. C9, one needs to be extremely careful with numerical precision. 
Round-off errors in the numerically produced power-law models as well as in the inertial-mode eigenfunctions may 
occult the fine cancellation and lead to artificially large coupling. 

TIDAL OVERLAP IN OTHER MODELS 

The derivation leading to equation (C5) assumes that the integrand /(O) is sufficiently smooth. What is 'sufficiently 
smooth' and in what situation does this assumption break down? It turns out that the break-down occurs for realistic 
planet models and that the tidal overlap is much larger than what one obtains for single power-law models. 

The smoothness assumption is violated if /(©) has a discreet jump inside the planet. Such a discontinuity is caused 
by the density discontinuity associated with a first-order phase transition region (e.g., gas-to- metallic hydrogen phase 
transition region at r/R ~ 0.80). Let the jump be A/ at O = Oo- It contributes a term, A/sin(nOo)/n ~ A//n, 
to the tidal overlap. Even if Af is small, this term may dominate for high order modes. Similar reasoning applies if 
/(©) exhibits a discontinuity at a higher order derivative, for instance, if the above mentioned phase transition is of 
second order in nature so that a discontinuity in the gradient of density exists. In this case, the contribution to the 
overlap integral is of order ~ A/'/n 2 . 

The smoothness assumption can also be violated if /(0) is infinitely continuous yet it (or one of its derivatives) 
has a sharp transition over a small region, namely, if this transition occurs over a width of A0 which encompasses 
only one node or less (An ~ n A0/7r/2 < 1). This can be caused by, e.g., a relatively sharp power-law index change 
inside the planet. As is discussed in §A, gas pressure inside Jupiter changes its nature from that of an ideal gas to 
that of strongly interacting molecules around r/R ~ 0.98. Here we observe a variation in the polytropic index over 
one pressure scale height, or over a thickness of Ar/R ~ 0.002. Within this narrow region, / (0) varies rapidly for an 

19 This is so because the functional value at one of the two boundaries (the equator) is not zero - sec Eq. (C5). 
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amount A/ , and f( 2 \<3) has a peak value of ~ A/ /AO. The overlap integral 



f ,2 c„ s („e) A e) d e - s """ 9)/(9) + e^glggj " 2 _ 
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tt/2 



cos(nO)/"(0)dO 



(Dl) 



can be dominated by the last term and yields A/' /n 2 if An <~ nAr/R <~ nAQ/ir/2 < 1, or n < 1/0.002 ~ 500. For 
n > 500, /(©) can be considered as sufficiently smooth and the analysis in Appendix C applies. 

We numerically confirm these conclusions by integrating /(©) cos(nO) using a range of density profiles. Here, we 
take /(6) = r 6 \J p sm f/ PR 2 /pdr/dQ, where p sur f is [1 — (r/R) 2 } 13 with the (3 value taken at the surface. This /(6) has 
the same asymptotic behavior as near both boundaries. 

We show that when a density discontinuity is superimposed to a single power-law model (dotted curves in Fig. D10), 
the overlap integral indeed scales as l/n. Also, if the model has a sharp (but continuous) transition in the value of (3 
over a radius of Ar (dashed and solid curves in Fig. D10), the integral scales as l/n 2 for n < Ar/R, while for higher 
n values, it behaves as is predicted by equation (C7). We have also studied integration results for two realistic Jupiter 
models taken from Guillot et al. (2004) (models B & D). Model D has a first-order phase transition (dotted curve in 
Fig. Dll) and so its overlap integral scales as l/n; while the same phase transition is considered to be second-order in 
model D, and the resulting discontinuity in density gradient (as well as the equation of state transition at r/R ~ 0.98, 
see Appendix A.l) causes the integral to scale as l/n 2 . 

When the density profile is not a single power-law (as is the case in this section) , we could not solve for inertial-mode 
eigenfunctions exactly. We could only obtain an approximate solution that is good to the second order in wavenumber 
(0(X 2 ), see Paper I). It is reasonable to suspect that the overlap results obtained by integrating such an approximate 
solution deviate from the true one. A definite answer to this suspicion will likely be provided by full numerical solution. 
However, we argue below that the deviation should be unimportant. 

The result of integrating a fast oscillation function, as is shown in this section and Appendix C, depends only on 
the boundary behavior and interior discontinuities in the envelope of such a function. It does not depend on the 
exact shape of the function in the interior. Our approximate solution to the inertial-modes is exact near the surface, 
and is sufficiently accurate near the center (where the WKB approximation works well). Moreover, when a density 
discontinuity (or discontinuity in density derivatives) is present inside the WKB region, as inertial-modes are insensitive 
to density structure, the solution is not expected to deviate qualitatively from the approximate solution that does not 
take this into account. 

In confirming the scalings derived in this section, we have only integrated the toy-model (/(OcosnO), instead of 
integrating inertial-mode eigen- functions (in Fig. C9 we integrate both). One can similarly argue that integrating an 
appropriately chosen /(©) is equivalent of integrating the real function. In fact, our toy model should produce results 
both qualitatively and quantitatively similar to that obtained using the actual eigenfunctions, one can almost make 
do without detailed knowledge of the latter.es such a density discontinuity. 

Lastly, independent of the radial profile, integration in the angular direction always introduces a factor of l/n 
cancellation. 
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Fig. D10. — The severity of cancellation in the overlap integral calculated using the toy model for three different density profiles are 
shown in the lower panel as a function of n (n even), while the top two panels show the corresponding /(©) (left) and df/d@ (right) as 
functions of f — r/R = 1 — cos0. We take /(©), the envelope of the cosine function to be /(©) = r e p 2 /p(p BUT f / p^/^dr /d& while the 
various density profiles are: a /3 = f .8 power-law model, overlaid with a 1% density jump at r/R = 0.8 (dots, exhibiting a 5-function in 
df/d&); a mock Jupiter model where the power-law index varies from 1 in the interior to 1.8 in the envelope, with the transition occurring 
at r/R = 0.98 (solid lines, having a jump in df /d&) and spanning a range of Ar/R ~ 0.002 (FWHM of the spike in d 2 f/d0 2 ); a similar 
model but with the transition occurring over a range of Ar ~ 0.02 (dashed curves, the one with smooth df/d@). Analytically, we expect 
scalings of 1/n, 1/n 2 switching to 1/n 3 - 6 when n > 500, and 1/n 2 switching to 1/n 3 - 6 when n > 50, for the three models, respectively. 
These scalings are marked here as the three dot-dashed lines. 
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Fig. Dll. — Same as Fig. D10 but with the density profile taken from two realistic Jupiter models: models B & D as in Guillot ct al. 
(2004). Model B (solid curves) is based on an interpolated equation of state with no core and no density discontinuity across the metallic 
hydrogen phase transition region at r/R ~ 0.8 - but the first derivative of density is discontinuous there (df/dQ jumps by ~ 50%). Overlap 
integral in model B is expected to suffer a cancellation with a 1/n 2 scaling (lower panel). The sharp transition in the equation of state 
around r/R ~ 0.98, with a FWHM for d 2 f/dQ 2 of Ar/R ~ 0.02, also contributes to this scaling. But this contribution falls off sharply 
for n 3> 1/0.02 ~ 50. Model D (dotted curves) has a 10M^ solid core, and is based on PPT equation of state with the phase transition 
being first-order, giving rise to a fractional density jump of ~ 20%. This is seen here as the jump in /(©) and the spike in df/dQ. Overlap 
integral in model D is dominated by the density jump and it scales roughly as 1/n, as expected. These results are insensitive to core sizes, 
since the r 6 scaling in /(©) near the center suppresses any influence from the inner boundary condition. Moreover, severity of cancellation 
calculated for actual inertial-mode eigenfunctions is expected to be one power of n steeper than those presented here, due to cancellation 
in the angular direction. 



